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THE EFFECT OF SUPPORT FLEXIBILITY AND DAMPING 
ON THE SYNCHRONOUS RESPONSE OF A SINGLE 
MASS FLEXIBLE ROTOR 

By: R. G. KIRK 1 

E. J. GUNTER 2 

ABSTRACT 

This paper deals with the dynamic unbalance response and transient 
motion of the single mass Jeffcott rotor in elastic bearings mounted on 
damped, flexible supports. 

A steady state analysis of the shaft and the bearing housing motion 
was made by assuming synchronous precession of the system. The conditions 
under which the support system would act as a dynamic vibration absorber 
at the rotor critical speed were studied and plots of the rotor and support 
amplitudes, phase angles, and forces transmitted were evaluated by the 
comp..ier and the performance curves were plotted by an automatic plotter 
unit. Curves are presented on the optimization of the support housing 
characteristics to attenuate the rotor synchronous unbalance response. 

The complete transient motion including rotor unbalance was examined 
by integrating the equations of motion numerically using a modified 4th 
order Runge-Kutta procedure and the resulting whirl orbits were plotted 
by an automatic plotter unit. The results of the transient analysis are 
discussed with regards to the design optimization procedure derived from 
the steady-state analysis. 
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INTRODUCTION 

The study of rotor dynamics has ir. recent years, become of increasing 
importance in the engineering design of power systems. With the increase 
in performance requirements of high-speed rotating machinery in various 
fields such as gas turbines, process equipment, auxiliary power machinery 
and space appl ’cations, the engineer is faced with the problem of designing 
a unit capable of smooth operation under various conditions of speed and 
load. 

In many of these applications the design operating speed is often 
well beyond the rotor first critical speed, and under these circumstances 
the problem of insuring that the turbomachine will perform with 3 stable 
low-level amplitude of vibration is often difficult to achieve. 

At the turn of the century H. H. Jeffcott (I ) developed the fundamen- 
tals of the dynamic response of the damped single mass unbalanced rotor 
on a massless elastic shaft mounted on rigid bearing supports. The Jeffcott 
analysis of the single mass model showed that operating speeds above the 
first critical speed were possible and that a low level of vibration would 
be attained once the rotor had exceeded the first critical speed. 

As various compressor and turbine manufacturers adapted the flexible 
rotor design concept in which the rotors were designed to operate above the 
first critical speed, various units developed severe operating difficulties 
which could not be explained by the elementary Jeffcott model. 

Under certain conditions of high speed operation above the first 
critical speed, such influences as internal rotor friction (2), hydro- 
dynamic bearing and seal forces* (3) and aerodynamic cross coupling (4) 
can lead to a destructive nonsynchronous precessive whirl motion being 
developed in the rotor system. 

B. L. Newkirk and Kimball (5), in their early investigations of self-excit- 
ed instability in compressors due to internal friction, were able to deter- 
mine experimentally that the introduction of a flexible support system 
could greatly extend the rotor stability threshold speed. D. M. Smith 
(6) in 1933 was the first to verify Newkirk’s findings theoretically by 
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expanding the Jeffcott model wirh internal damping to include a massless 
damped flexible support system. Recent investigators such as Gunter (7), 
Tondl (8), D i men t berg (9) and others (10) have shown that flexible damped 
supports may improve the stability characteristics of high speed ’-otors. 

The problem of bearing forces transmitted has been examined by various 
researchers, (II, 12, 13, 14). They have shown that a significant reduc- 
tion in the forces transmitted can be achieved by the proper design of 
the bearing support system. 

The present analysis was undertaken to determine the influence 
of flexible supports on the synchronous unbalance response of the single 
mass Jeffcott rotor, and to optimize the support system characteristics 
so as to minimize the rotor amplitude and forces transmitted over a 
given speed range. Der. Hartog (15) has shown that the tuned vibration 
absorber will greatly reduce the response of the forced vibrations of 
the two-mass system. The following analysis parallels this approach 
for the case of a single mass rotor excited by an unbalance load. 

This paper presents an analytic study of the tuned damper support 
system similar to that employed by Brock (16) and also presents a 
generalized study performed on the digital computer to obtain optimum 
support damping to produce the best response of the rotor over a wide 
speed range. It is well known that a damper support system can improve 
the vibration characteristics of a rotrting shaft and various investigators 
have considered the problem either from the standpoint of a continuous 
elastic system or as a series of lumped masses (17-23). 

Although the results presented in this paper apply specifically 
to the single mass Jeffcott model, the optimization procedure may be 
readily extended to more complex multi -mass rotor bearing systems by 
employing a finite element rotor digital computer program similar to 
the procedure presented by Lund in Ref. 24 or by using the procedure 
as outlined in the paper presented by Crook and Grantham (25) on the 
vibration analysis of turbine generators on damped flexible supports. 
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EQUATIONS OF MOTION 

Figure I Represents the single mass Jeffcott rotor mounted in damped 
elastic supports. In the Jeffcott model, the shaft is considered as a 
massless elastic member and the rotor mass is concentrated in a disc 
mounted at the center of the span. The shaft Is supported in linear 
bearings which are mounted in damped flexible supports. 


Neglecting rotor acceleration and the disc gyroscopics, the governing 
equations of motion for the rotor, bearings, and support system in complex 
notation reduce to the following (31) 


M 2 Z 2 + C g Z 2 + C.Z $ -|QZ 2 + (K s - iuC.)Z s = M 2 e u w 2 e ,u+ 


(I ) 


C.Z. - C.Z + K.Z. - (K - iu>C. )Z = 0 
b j is b j s is 


( 2 ) 


7i + CiZi + KiZi - C.Z - (K - iwC.)Z = 0 

1 11 1 * IS s IS 


(3) 


where 


I : 


Z = Z 2 - Z. - Zi = relative shaft deflection. 
s J 


If the internal damping C. and the aerodynamic cross coupling term Q 
are excluded from the above equations then the system will be stable (26). 

After the initial transient motion has damped out, it may be assumed 
that the system steady-state motion is circular synchronous precession. In 
this case the displacements are related to the velocity and acceleration 
vectors as follows: 


Z. = A.e 

i i 


iwt 


Z. - i«Z. 


( 4 ) 


Z. = itoZ. = - u> 2 Z. 

i i i 


* I I I ustrat i ons begin on page 60. 
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where A, is in general complex. 

The differential equations of motion may be reduced to a set of algebraic 
equations for the determination of the rotor steady-state motion. 


(K s " 

M 2 uj 2 + iC w)A 2 - K A. - K Ai = M 2 e to 2 
c s *■ SJ S * 2 u 

(5) 

-k s a 2 

+ (K. + K + iwC. )A. + K A, =0 
b s b j s 1 

(6) 

-k s a 2 

+ K g A. + (Kj + K s - Mjo) 2 t itoCj )Aj = 0 

(7) 


ROTOR AMPLIFICATION FACTOR 


Consider the steady-state orbit of the flexible rotor on rigid 
supports. The rotor amplitude is a function of both the rotor and bearing 
stiffness and damping characteri sties. Assuming Aj is zero, the relative 
journal bearing complex amplitude from Eq. 6 is given by 


K (K + K w - io)C, ) 

A. 3 6 6_A2 

1 (K *- K, >2 + (uCJ 2 
so b 

Solving Eq. 5 for the rotor amplitude yields 

( K2 “ I CjjC 2 5 

A2 = 

“ U (K 2 - M 2 u 2 ) 2 + (a)C 2 ) 2 


( 8 ) 


(9) 


where 


Ko = 


K, K (K + K, ) + K ( u>C , ) 2 
b s s b s b 

(K + K.) 2 + ( 10 C .) 2 
s b b 


C 2 = 


K 2 C. 
s b 


(K. + K ) 2 + (wC, ) 2 
b s b 


+ C 


The rotor displacement vector Z 2 may be expressed in terms of the 
absolute displacement R 2 and the "base angle t as follows 


Z 2 


~ R 2 e 


i (wt - <f> ) 


( 10 ) 
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( 


where 


R, = 


l^e^w 2 


~\JTK 2 ~ M 2 to 2 )^ + ((0C2 
U)C 2 


<j> = tan -1 


K? - M 2 w 2 


The above results are similar to the rotor amplitude and phase angle 
results for the single mass flexible rotor on rigid supports as shown by 
Thomson (27). 

The rotor undamped, or natural critical speed is given by 



For the case of a lightly damped rotor system on rigid supports the 
maximum rotor amplitude will occur at approximately the rotor critical 
speed and the dimensionless rotor amplitude or amplification factor at 
the critical speed is given by 

R 2 k 2 

A - — | = — rr- (12) 

e ‘ to=a) to C 9 
u c c *■ 


Example I 

Consider a 97 lb. disc centered on a uniform massless elastic shaft 

K k 

as shown in Fig. (I). Assume that the bearing stiffness ^ is 500,000 
lb/ in and that the effective shaft stiffness K g at the disc station is 
333,000 Ib/in. Assuming light damping, the total stiffness K 2 is given by 


«2 = 


K K. 
s b 


K + K. 
s b 


I x 0.333 x IQ 12 
(I + 0.333)1 0 6 


= 250,000 Ib/in 


The rotor cri + ical speed is 


0) 


c 


_-\ /** _~\ / 250,000 
V M 2 ~ V 0.25 


1,000 rad/sec 
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or N = 9,550 RPM. 
c 

If the rotor damping C g is assumed to be 15 lb - sec/ in and the bearing 
damping coefficient C^/2 is 80 lb - sec/ in. fhen the effective system damp- 
ing coefficient C 2 is approx ima+e I y given by 


C 2 =C - 15 + (0 ; 333ii JL j£i iL l60 = 25 . sec/in 

" (K. + K ) 2 (I.333) 2 x I0 12 

b s 

The amplification factor at the rotor critical speed is given by 
K 2 


CR 


A = 


250 ,000 . = 10 0 

r\r\r\ oc I *\J 


a> c C 2 1,000 x 25 


The amplification factor of 10 represents a very lightly damped rotor 

system and indicates that the rotor amplitude at the critical speed will be 

10 times the rotor unbalance eccentricity e . 

1 u 


ROTOR RESPONSE ON DAMPED FLEXIBLE SUPPORTS 


Solution of Eg. 6 for the case of synchronous precession for the 


snaft relative deflection Z yields 

r K, (K, + K ) + (wC. ) 2 + iwC.K 
x I b b s b bs 

/ - (Z 2 - Zj) 

(K. + K ) 2 + (ojC, ) 2 
L b s b 

Hence, in terms of the general coefficients C 2 and K 2 


( 13 ) 


Z s = (Z 2 - Zj) 


K 2 + iu)(C 2 “ C ) 


K 


(14) 


The simultaneous equations for the absolute shaft ana support housing 
motion reduce to the following 


TKj + Ko * + 1 (x) ( C 2 


C flAj + C-iu)(Cp * C ) * Z]A 2 = 0 


(15) 
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C“K 2 - i u (c 2 - C S )?|A! + CK 2 - M 2 U 2 + iuC 2 ]A 2 ' (16) 


If the damping terms are neglected, then the natural frequencies 
of the system may be determined by the expansion of the determinant of 
coef f icien'i s. The resulting frequency equation may be exp^ssed ?s follows: 


= Vl + i 2r !S± ~\f ( 


JLJS + 1)2 . is 

2M 2 M 


(17) 


where 

u> 

C 

Figure 2 re -esents ihe dimensionless critical speeds vs. the dimen- 
sionless support stiffness factor K for various values of support to rotor 
mass ratios. Mote thaT the incorporation of the Hexible support with the 
rotor bearing system causes two critical speeds to occur; one which is 
higher and one which is lower tnan the original rotor critical on rigid 
supports. 

To solve for the complex support and rotor amplitudes Aj and A 2 , Eq. 

1 5 and 1 6 may be expressed as follows: 



[a.. + ib.jA. = F.; j = I, 2; ».= I, 2 (18) 

ij ij j ; J 

Multiplying Eq. 18 by the complex inverse matrix of coefficients and ex- 
panding yields 


F j ai 2 + ib 12 

F 2 a-> 2 + > b 22 


(19) 


Where 


A = d + id. 

r I 

d^_ = ( K 2 — M 2 oj2 ) (Kj - Mjuj2) — K 2 M 2 (u 2 — CjCjiU) 2 — to 2 Cg(C 2 - C^) 
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d. = C 1 u(K 2 - M 2 u 2 ) + C 2 “(K 1 - M lU 2 - M 2 u 2 ) + C s w(K 2 + M 2 u 2 ) 


Expanding Eq. 19 

Fia 22 - F 2 aj 2 + i(Fl&22 “ ^2^12^ 


Ai = 


d p + id. 


(20) 


In this case only an external unbalance excitation force F 2 is acting 
on the shaft and no external exciting force F i is assumed to be present on 
the support syster. For example, an excitation force Fj may be transmitted 
to the rotor system through the support structure by vibrations of auxiliary 
or adjacent equipment. 


Ai = - 


F 2 Cai 2 c lr + ^12^ j + Hbi 2 <l r - a^dj)!] 


d 2 + d. 2 

r i 


( 21 ) 


Assume is of the form 


A l A ir - i . 


( 22 ) 


The complex support amplitude Zj after some complex algebraic manipulation 
is given by 


2, = A,e'“ + = R,e i< “ + ‘ 8,1 


( 23 ) 


where 


d. 


R l = V* A ir 2 + A li 2 * 61 = +an_ll T* ) 


If the shaft damping coefficient is considered small in comparison 
to the effective damping coefficient C 2 than the system displacements and 
phase angles are given as follows 


/K 2 2 + ( ii)C ? ) 2 

Ri = M ? e uj 2 ' 

U »d r 2 *1. 2 


( 24 ) 


and the phase angle of the support motion relative +o the rotating unbal 
,.'d is giv. n p / 
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1 


I - 


$1 = tan 


-l 


l<2d . ~ wC2<l r 


K 2 d r t u>C 2 d. 


( 2 ?) 


Since the complex rotor support motion Zj is given by 

i (ut - 8i ) 

Zi = Xi + i Y i = Rje 

Then for example, the horizontal and vertical components of the support 
motion are given by 


| / K 2 2 + (uC 2 ) 2 $cos(wt - Si ) ^ 

" - (d 2 -\ / 

r i 




i =M 2 e. 


2 \ Jl \ K 

y d 2 + d 2 |si n(u»t - Si)) 


(26) 


In a similar fashion, the complex rotor amplitude Z 2 is given by 


ail + 5b li 

Z, = M 2 e u d + i d . e 
r i 


iut 


(27) 


After some manipulation, Eq. 27 reduces to the following 


Z 2 = M2 e u “ 2 ' 


/(Kj + K 2 ~ Mjai 2 ) 2 + FTCi + C 2 )u) 2 i (ut - 3 2 ) 
e (28) 


d 2 + d. 2 
r 1 


where 


/ («! + K 2 - M lW 2 )d. - (C l + C 2 )ud 

e 2 = tan" 1 j ! - 

^ (Kj + K 2 - Mjio 2 )d r + (Ci + C 2 )u>d . 

The relative journal displacement is given by 


(29) 


2 j ■ Z ’ - z ‘ - Z S 


(30) 


Where the relative shaft deflection is 
(Z 2 - Zj) 


Z = 
s 


K 


* K 2 “ iu)C 2 ] 


(31 ) 


Solving for the journal displacement 


1 


• j 



t 


* 

l 


— * ■ -r 


4 . 


f 

* I 

I 

I 


" 1 


\ (cot - 6 . ) 

Z. = R. e J 

J J 


where , 


R. = M 2 e u> 2 - 
J u 


‘ (Ki - MjX) 2 + (uCj ) 2 " 

v 


f K s- K f + ] 

juC 2 ) 

2“ 

d 2 + d . 2 
r i J 

A 


I K s ) 1 

i K si 

. 


(32) 


and the phase angle Bj between the journal amplitude and rotating un- 
balance force is given by 


6 - 

J 


tan -1 


/ (Kj - Mjw 2 )d . - uC i d r \ 
\ “ Mjo 2 )d. + diCid.f 


+ tan -1 



(33) 
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FORCES TRANSMITTED 

The magnitude of the resultant forces transmitted through the bearings 
and the support are of considerable interest to the designer from a stand- 
point of bearing life and system isolation. It is desirable to minimize 
the forces transmitted through the supporting structure and foundation so 
that other machines or piping systems are not excited. The magnitude of the 
force transmitted through the bearings ;s given by 

F. = R. -v/k" 2 + (wC. ) 2 (34) 

o j v b b 

and the force transmitted through the support system is given by 

F! = Rj VKi 2 + (wCj) 2 (35) 

An indication of the effectiveness of the support system ’ jftenuating 
the forces transmitted to the foundation is the support dv..imic transmissi- 
bility factor TRD which will be defined as the ratio of the magnitude of the 
transmitted support force to the rotating unbalance load. If the dynamic 
transmissibi I ity is less than I, then the support system possesses good 
attenuation characteristics. Analysis has shown that if the support 
housing impedance characteri sties, which are determined by the housing mass, 

10 





stiff' 3ss and damping, are mismatched to the rotor-bearing system then 
i nde- certain speed conditions the dynamic transmissibi I ity may exceed I. 


The dynamic transmissibi I ity for the support is defined as 


i'RD 


F j / ( K 2 2 + (uC 2 ) 2 )(K 1 2 + 

M 2 e u» 2 V d 2 + d. 2 


(C^) 2 ) 


(36) 


f it is assumed that the rotor is operating well above any of the 
sysre i critical speeds then the dynamic transmissibi I ity is approximately 
given by 


TRD 




(K 2 2 + (wC 2 ) 2 )(K 1 2 + (C^) 2 ) 
Mj 2 M 2 2 


(37) 


The above expression leads to the well known conclusion tr.at to 
mir imize the forces transmitted through the suppor J for supercrit ical 
sp< ed operation in the Jeffcott model, the support damping should be 
ze n and the support stiffness should be as light as possible (28). 

This is i highly ...'desirable design practice for several reasons since 
large o!or amplitudes ana forces 1 ra.ismitted may be encountered when pass- 
ing tnrotgh the rotor critical speeds, and also the rotor system would be 
extreme’y shock sensitive and particularly susceptible to self-excit.d 
whirl instability under such conditions. 

A compromise support damping coefficient should oe selected to either 
minimize the rctcr amplitudes or the forces transmitted over the operating 
'ed range a also be sufficient to insure adequate rotor stability. 


ANALYSIS OF SYSTEM UNBALANCE RESPONSE - TUNED SYSTEM 


Figure 3 renr^sents a computer generated plot of tne dimensionless 
rotor relative amp'itude versus the dimensionless rotor speed for the 
case of =-i = I . This relative rotor amplitude is equivalent to the 
motion monitored by a proximity probe mounted in the casing measuring 
the rotor motion at the center span. This system represents a tuned 
condition in which the support stiffness ratio K is equal to the support 


1 1 


- *•* 


i 
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mass ratio M. With no support damping in the system, the tuned support 

will cause the relative rotor amplitude to be zero at a speed correspond- 

I 

5 ing to the rotor critical speed with rigid supports. The introduction 

of support mass and flexibility has caused two critical speeds to appear 
in the system; one above and one below the rigid support rotor critical. 

Note that when the support damping is relatively low the amplitudes at 
the two criticals becomes extremely high. 

As the dimensionless support damping ratio C increases fron. 0.01 to 
10 The rotor amplitudes at the system critical speeds decrease while the 
amplitude increases at a speed corresponding to the rigid support critical 
speed (w/o) c = I). Note that in this case the damping value of 10 appears 
to be close to an optimum value for the minimization of the resonance 
amplitudes. If the support damping is further increased from 10 to 50, 

Fig. 3 indicates that there will be only one critical speed present in 

the system which will correspond to the rigid support critical. Although 

the damping of C = 50 is over 5 times the optimum value, the maximum 
amplitude is only 1/3 the rigid support value of 10. As the support damp- 
ing approaches infinity, the rotor amplitude will asymptotically approach 
10 . 

Figure 4 represents the absolute dimensionless rotor motion for various 
values of support damping ratio and is similar to Fig. 3. It should be 
noted that the damping coefficient of 10 also appears to be close to the 
optimum damping for the absolute motion as well as the relative motion. 

It is of interest to note that The various damping lines all intersect 
at a common point P in the plot of absolute as well as relative rotor motion. 
If the rotor amplification factor A is 100 (implying light rotor damping) 
then there will be two common points of intersection P and 0 on The response 
plots (see Fig. 10) similar to that shown by Den Hartog for the damped 
vibration absorber (15). The intersection points P and Q will occur at speeds 
respectively below and above the rigid support critical speed. The rotor 
amplitude may be minimized for the case of the absolute rotor motion by 
f selecting the damning such that the slope of the response curve is zero 

at point P, and zero at point Q to minimize the rotor relative mofion. 
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Figure 5 represents the phase angle beTween the rotating unbalance 
vector and the absolute rotor displacement vector for various damping 
coefficients. The phase angle for the single mass rotor on rigid 
supports (Jeffcott model) increases with speed from 0 to 90 degrees 
at the critical speed and asymptotically approaches 180 degrees as the 
rotor speed greatly exceeds the critical speed. The phase angles of the 
rotor on damped flexible supports has a considerably different behavior 
from that of the rigid support rotor. For light values of support damping 
(C = 0.01), the phase angle increases rapidly to 180° as the system passes 
through the first critical speed and drops to almost 60° as it passes 
through the second critical speed. As the speed greatly exceeds tne 
highest critical speed, the phase angle again approaches 180°. The phase 
angle of 180° indicates that the rotor mass center lies along the rotor 
spin axis. As the support damping coefficient is increased beyond 5 for 
the case of the tuned system, the reduction in phase angle above the first 
critical speed is suppressed. This phenomena of phase angle reversal 
above the first critical speed has been observed experimentally (30). 

Figure 6 represents the support amplitude versus speed for various 
damping values and indicates that with very light support damping there 
will be large support resonances. As the damping is increased beyond 
C = 10 the resonances are suppressed and the amplitude is only slightly 
greater than I. For C = 50 there is only a small peax observed in the 
support system which occurs at a speed corresponding to the rigid support 
critical speed. The addition of high damping (C > 50) freezes the support 
and limits its motion drastically. 

Figure 7 represents the support phase angles versus speed ratio for 
various values of support damping. The phase angle for light damping 
(C = 0.01) is zero at low speeds and goes to 180 degrees as it passes 
through the first critical and then shifts to 330° upon passing through 
the second critical speed. If the rotor damping is light (A = 100) the 
support phase angle will approach 360° after passing through the second 
critical speed. Note that the various damping lines intersect at three 
points. The first node point represents the first system critical speed, 
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the second node point represents the rigid Support critical speed and the 

» 

third node point represents the second critical speed on flexible supports. 
In the discussion of the single mass flexible rotor presented in vibration 
texts (27) the phase change is only shown from zero to 180 degrees. In 
more complex systems with flexible supports, the ph-’se change may vary 
between 0 and 360 degrees. For example in multimass systems the authors 
have observed phase changes of n times 180 degrees where n represents the 
number of system critical speeds. The measurement of rotor and support 
phase angles have been neglected and limited data has been reported in the 
literature. This is an extremely useful variable which when xorporated 
with displacement measurements can be used in balancing flexible rotors 
or impedance calculations of the support system. 

Figure 8 represents the dimensionless bearing forces transmitted for 
the tuned system. The dimensionless force transmitted is obtained by 
dividing by the transmitted force corresponding to the value at the crit- 
ical speed of the original rotor on rigid supports. Because of the light 
shaft damping the force transmitted curves are similar in appearance to 
the displacement curves. Note that for the support damping coefficient of 
C = 10 the forces transmitted to the bearings at the rigid support critical 
are only 10 percent of the value transmitted for the rotor bearing system 
on rigid supports. 

Figure 9 represents the force transmitted through the bearing supports 
to the foundation or base for various values of supporting damping. With 
a very lightly damped support system, (C = 0.01) the support amplitude and 
force transmitted will be particularly high at the first critical speed 
where the bearing and support motions are in phase. At the second critical 
speed, the support amplitude is lower than the amplitude attained at the 
tfrst critical speed. This is because the bearings and support motions 
are out of phase which enables the bearing damping to help attenuate the 
support motion. It is of interest to note from Fig. 8, for the tuned 
rotor system, the bearing force transmitted at (u/u> c ) = I with an jidamped 
support system is zero. Figure 9 shows that the correspond i ng force trans- 
mitted through the support system at w/u) c = I has been reduced to only 10^ 
of the rigid support value. 
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The force transmitted for an undamped support system at a speed 
ratio of four is approximately 10? of + he rigid support value. This con- 
dition would be desireable if it were possible to accelerate through the 
criticals, thereby avoiding the large steady-state amplitudes and forces 
developed. 

The near optimum damping of 10 increases the support forces transmitted 
in the supercritical speed region to 30? of the rigid support value and 
the overdamped support system (C = 50) has increased to ?»H y 80?. Hence, 
the support damping introduced to suppress the system resonances will cause 
the forces transmitted to increase in the supercritical speed region. 

If the system is designed to operate over the entire speed range shown, 
then the near optimum value of damping (i.e., 0 = 10) for suppressing the 
rotor absolute amplitude also produces the most desirable attentuation of 
forces to the system support structure. 

OPTIMUM DAMPING FOR TUNED SYSTEM 

From the observation of the computer generated displacement and force 
transmitted plots it is apparent that there exists an optimum damping to 
either minimize the rotor amplitudes or the forces transmitted over the 
entire speed range. 

For example to minimize the absolute rotor motion as shown in Fig. 4 
or the relative rotor motion shown in Fig. 3, the method of (16) may be 
used in which the damping is chosen so that the slope of the amplitude 
curve is zero at points P and Q respectively. In the tuned system where 
K/M = I for light rotor damping (A = 100), the rotor amplitudes at points 
P and 0 are independent of the support damping as shown in Fig. 10 and 
can be shown to be equal to 

X 2 = x 2/e u l p> p = /l' + 2M (38) 

Therefore with the tuned system illustrated with a mass ratio of M = I, 
the maximum amplitude at P or Q will be 1.732 times the rotor unbalance 
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eccentricity. The optimum damping may be selected so that the tangent to 
the amplitude curve at either point P or Q has a zero slope. By selecting 
the optimum damping in this fashion it is seen that the maximum amp I itude 
in the system will not exceed the value given by Eq. 38. Thus it is readily 
apparent that to minimize the rotor response over a given speed range, the 
support mass should be kept as light as possible. 

After considerable algebraic manipulation (28) the optimum damping 
coefficient for both points P and Q is given by the following expression 

52 _ 4M 3 ij > 3 - 3M 2 (4M + 3 )» 2 + M(I2M 2 + I3M + 8 )^ - M(l + 2M ) 2 (3g) 

-1 2M<Ji 2 + 8(1 + 2M)ip - (I + 2 M) 
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where. 


“c 

€ = Cj/C c = C x /C 2 x I/2A = Cj 

1 i> = fy 2 or £ 2 2 2 depending on whether the value calculated is for point 
P or Q respectively. 


and 


£2 


^1 + 2M 

I + /TTlM 


„ ? _ A + 2M 

“2 ~ , 

A + 2M - I 

For example, when M = I and for the first node, P: 


* » fli- « — = 0.634 

I + A5 


and 


? 2 = 0.447 


Hence 
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— = 0.688 for point P 

c opt 

In a simi lar fashion 
Ci 

■x-l = 0.559 for point Q 

c opt 

Example 2 

As an example of the application of the tuned support design criteria 
consider the rotor of Example I mounted in flexibly supported bearing housings 
which weigh 48.5 lbs and have a stiffness of 125,000 Ib/in. The total 
support weight Wj and stiffness Kj is given by 

Wi = 2 x 48.5 = 97 lb 

Kj = 2 x 125,000 = 250,000 Ib/in 

Hence, 

M = Mi/M 2 = I .0 
K = Ki/Ko = !*0 


The critical damping coefficient C is given by 

c 


c = 1 Z . SOO^OOp \y±i. - 500 lb-sec/in 
c w 1,000 rad/sec 


Thus the support damping coefficients required to make the slope of 
the rotor amplitude curve zero at points P and Q are respectively given 
as fol lows 


Ci I = 0.688 x C = 344 lb-sec/in 
1 Ip c 

Ci I = 0.559 x C = 279.5 lb-sec/in 

1 'q c 
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These calculations are valid only for the case of zero damping on 
the rotor and in the bearings (i.e., A = ») and only for the tuned system 
(i.e., K = M). For a more realistic solution, a value of A = 10 was chosen 
and numerous cases were then programmed on a digital computer to arrive at 
a value of optimum amplitude and required damping. This approach is dis- 
cussed in the next section of this paper but the results for the tuned 
system are very nearly the same as the results arrived at analytically for 
the case of A = « and are presented in Fig. il. 

The results shown in Fig. il are approximately correct for systems 
having moderate to light damping on the rotor (i.e., 10 £ A < »). Note 
that the smaller the mass ratio M, the lower will be the peak response 
and also the lower will be the required support damping. For example, if 
the mass ratio is 0.1, then the maximum dimensionless amplitude will be only 

I. l and the required damping ratio will be 5 as compared to a value of 13.6 
for an M ratio of I. Figure 12 is a response plot for the tuned system 

K = M = 0. I which ill ustrates the val idity of the results plotted in Fig. 

II. The response curve for a damping ratio of 5 passes almost horizontal 
through the node point and has the low amplitude ratio as indicated by 
Fig. I I . 

Example 3 

Consider a rotor system similar to Example 2 in which the rotor rigid 
suoport amplification factor A = 10. 

For a tuned support system the dimensionless support damping coefficient 
is obtained from Fig. II for M = I as follows 

C = C l /Z 2 = 13.6 

where C 2 is given as 25 lb-sec/ in (Example I). 

Therefore, 

Ci = 13.6 x C 2 = 340 lb-sec/ in 
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Note that this value is approximately the same as the value given in Example 
2 for the required damping at point P corresonding to A . 

This indicates that each support must have 170 I b-sec/ i n. damp ing to 
achieve the optimum response of about 1.7 times the unbalance level of the 
rotor. 

Next consider a tuned support with a mass u.id stiffness ratio of 0.10 
(see Fig, 12). Corresponding support weight and stiffness are given as 
fol lows 

Wi = 9.7 lb/ in 

Kj = 25,000 lb/ in 

The required damping is thus found from Fig. || to be 
C = 5.0 


or 


Cj = 5 x 25 = 125 I b-sec/ in 

Thus only 62.5 I b-sec/ i n. damp ing per support is required to obtain an 
optimum response of l.l times the unbalance level of the rotor. This value 
of l.l is in comparison to a maximum response of 10 times the unbalance 
level for the rigidly mounted rotor-bear i ng system. 

OPTIMIZATION OF SUPPORT DAMPING FOR OFF-TUNED CONDITIONS 

In general it is not possible or necessarily desirable to have a 
tuned support system. The support to rotor mass ratio is usually dic- 
tated by design considerations and can be varied only within certain 
ranges. Figure I I shows that for best reduction of rotor amp I itude, the 
support mass should remain as light as possible. However, it will be 
shown that even with high mass ratio support systems the rotor ampl itudes 
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can be attenuated by a factor of 5 by proper selection of the stiffness 
and damping coefficients. 

To evaluate the optimum damping for off-tuned conditions the computer 
program was run for various support mass and stiffness ratios and each of 
these for various damping coefficients. For example, Fig. 13 represents 
the amplitudes at the rotor first and second critical speed for various 
mass ratios with a dimensionless stiffness ratio of K - 0.01 as the mass ratio 
and damping are varied. The solid lines represent the amplitude at the 
second critical speed and the dotted lines represent the ampliiude at the 
first critical speed. With moderate support damping ratios it is observed 
that as the mass ratio increases the amplitude at the first critical reduces 
while the amplitude at the second critical increases. The optirm damping 
was selected as the intersection of the amplitudes at the first and second 
critical for a particular value of damping. For example in Fig. 14 for K = 1.0, 
the lowest optimum amplitude poin+ on the p I o r 's given by a damp'ng ratio of 
10 and produces an amplitude ratio of about 1.5. Several plots similar 
to Fig. 13 were produced and the results were then crossp lotted to obtain 
plots of amplitude versus damping ratio such as Fig. 15 for K = 5.0. 

Figure 14 represents the maximum rotor amplitude vs. support damping 
ratio for various values of dimensionless support stiffness for a rotor 
bearing system with a low support mass ratio of 0.01. Figure 16 shows 
that for this particular case, the lowest amplitude is achieved by a low 
support stiffness ratio of K = 0.01 which is of the same order as the mass 
ratio. With this low support stiffness, there is a wide range of support 
damping (i.e. C = I 6) that can be used to achieve the low 'eve I of 

rotor response. 

Thus, under proper design conditions the suppori damping may be 
allowed to vary by a considerable amount without impairing the rotor 
performance. As the support stiffness ratio increases, The maximum rotor 
amplitude response also increases and the required support damping must 
be larger. For example, if the support stiffness ratio increases from 
0.01 to 2.0, the optimum damping required increases by a factor of 15 from 
approximately 2 to 30. 
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Figures 17 and 18 represent the maximum rotor amplitude for mass ratios 
of 0,10 and 0.715. Note also that for high stiffness ratio support systems, 
the permissible range of the support damping coefficient is very narrow, 
and that either a reduction or an increase o r damping beyond the optimum 
value will result in a rapid gain in rotor response. 

it is also of interest to note that if a high support stiffness (K = 2) 
is used in conjunction with a low value of support damping (C < 2) then the 
rotor response will be worse than the original rotor response on rigid 
supports vA = i0). 

Figure 19 represents The maximum rotor response vs. support damping 
for a high support mass ratio system (M = 2). It is obvious from tne 
comparison of Figs. !6 and 19 that the high raass ratio support sysTem is 
less desirable. The minimum rotor amplitude that can be achieved is x 2 ^ e u ~ 2 
with a tuned support where K = M = 2 and a support damping coefficient of 
C - 20. (Also see FJq. II on the tuned system.) As the support stiffness 
ratio is reduced, the rotor response curve increases in the optimum damping 
region. 

If it is not possible to incorporate a high value of support damping 
into the system (C = 20), then the rotor amplitude can still be reduced 
to 40? o £ rhe original rotor response by a low support damping value of 
C = I and a reduced support stiffness ratio of K = 0.7. For low values of 
support damping, if the support stiffness increases beyond K = 0.7, the 
rotor response rapidly increases. 

A series of plots similar to Figs. 16 - 19 were produced for various 
mass ratios in order to determine the optimum rotor response for off-tuned 
support conditions. Figure 20 represents the rotor maximum amplitude vs. 
the support mass ratio for various values of support STiffness with optimum 
damping. 

For the case of A = 10, Fig. 20 illustrates that the lowest amplitude 
can be achieved with a low mass ratio support system. With a high mass ratio 
support system such as M = 5, the r o+or amplitude X 2 can be reduced from 10 
to 2.8 by moans of a tuned support stiffness of K - 5.0 and optimum damping. 
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Note that as the support stiffness becomes very light, the maximum rotor 
amplitude increases to 7.5. 


At a low value of support mass (M = 0.1), the rotor amplitude increases 
as the support sTiffness increases. The optimum damping required with the 
tuned supporT is given by the following approximate relationship 


1 . 37 x «2 0 .437 

c j x M 


u) 


(40) 


Figure 21 represents the rotor maximum amplitude vs. stiffness ratio 
K for various values of mass ratio M and damping val jes C. This figure 
illustrates that it is possible to operate with off-tuned conditions and 
still maintain a low level of vibration. It is seen that the lignt damping 
value C = 0.10 will produce the highest amplitude over the range of 
stiffness plotted for K = 0. I to 10. For K values less than 1.5, the damping 
ratio should be less than 10, while for high stiffness supports where K > 2, 
the damping value C should be > 20 for maximum attenuation. Note ihat for 
low stiffness supports (K < 0.2) the value of C = 20 represents an over- 
damped support system causing the amplitude and transmitted forces to be 
greater than the optimum value. 

Figure 22 represents the rotor maximum amplitude vs. support damping 
ratio C for various values of support stiffness K with optimum support 
mass M. Figure 22 shows that the lowest amplitude level can be achieved 
with a low support stiffness (K = 0.01) and the support damping may vary 
from C = 0.5 to 10 while maintaining a low level response. The figure also 
illustrates that as the support stiffness is i ncreased, the amplitude will 
also increase for a given value of damping. It is also clearly seen that 
as the stiffness value is increased, a larger support damping value is re- 
quired to produce a low vibration amplitude with optimum support mass ratio. 

Figure 23 represents optimum damping and mass ratios for various values 
of stiffness ratio. Figure 23 shows that, as the mass ratio increases, the 
required stiffness ratio increases for a given value of damping. It is of 
interest to note that for K values between 0.2 and 2.0, for a given M value, 
there can be two values of optimum damping, a low value of C below !0 and a 
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high damping va'ue > 10. Although high damping values may result in low 
rotor amplitudes, the bearing forces transmitted through the support will 
be much higher. Therefore extremely large values of support damping 
should be avoided. 
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TRANSIENT ANALYSIS 


The previous discussion has been concerned only with the steady-state 
response of the ro v or due to unbalance and has not considered the rotor 
initial transient notion. As discussed previously, the damped flexible 
support system is important, not only from the standpoint of reduction of 
synchronous unbalance response, but also in the control of self excited 
vibrations such as caused by internal friction, aerodynamic excitation, etc. 
Therefore to investigate the general rotor motion and also to provide a 
check on the steady-state analysis, the rotor equations of morion were 
integratea forward in time on the digital computer using a modified 4th 
order Runge-Kutta integration procedure. This procedure is of importance 
particularly if the analysis is extended from a linear bearing or support 
system to include a nonlinear hydrodynamic damper bearing as presented in 
Ref. 13. 

The dimensionless rotor and support transient orbits were automatically 
computer plotted with the following dimensionless parameters 

X = x/e u . Y = y/e u 

Figure 24 represents the initial transient orbit of a 96.6 lb rotor 
of Example I with a highly damped supoort (C = 43) for the first 12 
cycles of shaft motion. The support mass ratio and the support stiffness 
ratio are both approximately the same (0.10) which represents a tuned 
system. Because of the excessive support damping, the maximum force 
transmitted to the support is 2.16 times the unbalance force while the 
force transmitted to the bearings is reduced by about 40%. The magnifica- 
tions of the force to the support would be highly undesirable for applica- 
tions such as aircraft jet engines. For example, various investigators 
have observed that such a situation occurs with the hydrodynamic squeeze 
film bearing when oporating at excessive eccentricity ratios (29). 

Figure 25 represents the bearing absolute and relative motion corresDond- 
ing to the case as shown in Fig. 24. The solid line represents the bearing 
absolute motion wni le the dashed line represents the bearing relative 
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motion. Since the support damping is 10 times the bearing damping, the 
initial absolute bearing motion is not much larger than the bearing rela- 
tive motion. Note that the timing marks on the orbit appear in the nega- 
tive x direction This indicates that the bearing motion is 180° out of 
phase with the rotating unbalance load. 

figure 26 illustrates the support housing motion. Because of ex- 
cessive support damping, the initial support transient motion is quite 
small and is less than the unbalance eccentricity. The pnase angle between 
the support motion and the rotating unbalance load is approximately 220°. 

Figure 27 represents the transient orbit for the same rotor system 
except that the support damping has been reduced by a factor of 100 from 
Cj = 1,000 Ib-sec/in. to 10 Ib-sec/in. In this case, the maximum force 
transmitted through the support is less than 9 % of the rotating unbalance 
force and the bearing force transmitted is 16?. This orbit is analogous 
to a suddenly applied unbalance such as a blade loss in an engine. Although 
the forces transmitted have been greatly reduced with the low ^riffness 
and damping support system, the rotor has developed a large in : tial transient 
motion of over 10 times the unbalance eccentricity and this transient 
motion is not readily damped out. 

Figure 28 represents the absolute and relative bearing moTion with 
the low support damping of Cj = 10 Ib-sec/in. The absolute bearing 
initial transient motion is extremely large while the relative motion is 
well behaved. Note that tne bearing relative phase angle has shifted from 
180° for the highly damped case to 50° for the case with light support 
housing damping. 

Figure 29 represents the support housing motion corresponding to the 
system with light support damping. A comparison of the absolute support 
motion and the absolute rotor motion indicates that the two are similar. 

This implies that the initial transient motion of the rotor is c;;« primarily 
to the large deflections in the support system. 
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In Fig. 30, the rotor transient motion is depicted with an optimum 
damping coefficient of C = 5.5 for minimum rotor response as determined 
from the steady-state analysis. The transient response is rapidly suppressed 
after seven cycles of shaft motion to produce a smal I stable synchronous 
orbit. The transmitted forces to the bearings and support are nearly 
balanced to achieve approximately a 75? attenuation of the unbalance load. 
Figure 30 shows that with the optimum damping as determined by the steady- 
state analysis, the initial transient rotor motion will be 5 times the rotor 
unbalance eccentricity. 

Figure 31 represents the bearing m + ?on. After approximately 6 cycles, 
the initial transient motion is damped out. Figure 31 indicates that the 
absolute bearing motion is equal tc the rotor unbalance eccentricity after 
the transient has died, and is 180° out of pnase with the rotating in- 
balance load. The relative bearing motion is approximately 60? of the 
unbalanced eccentricity and lags the rotating unbalance by about 120°. 

Figure 32 represents the support motion anc is also similar to the rotor 
motion as shown in Fig. 30. 
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SUMMARY AND CONCLUSIONS 

The equations of motion for a single mass rotor-bearing system on 
damped flexible supports have been derived and studied considering both 
a steady-state and transient type analysis. Design charts for both tuned 
and off-tuned support conditions have been presented. 

The analysis may be summarized by the following general statements. 

1. The critical speed response of the single mass Jeffcott model 
rotor may be completely eliminated by means of a low mass ratio flexible 
support with optimum damping. In this case the rotor steady-state ampli- 
tude of motion over the entire speed range will only be slightly mere than 
the rotor unbalance eccentricity. 

2. The support mass ratio should be kept as light as possible to 
achieve minimum rotor amplitude. 

3. Tne rotor amplitude may be considerably attenuated even for 
high mass ratio support systems by tuning the support stiffness such that 
K = M and incorporating optimum damping for the tuned conditions. 

4. With a low mass ratio support system, the required value of 
optimum damping is not critical and can vary by a factor of 10 without 
appreciably effecting rotor performance. As the mass ratio increases, the 
required value of optimum damping increase-:'' rapidly and the permissible 
range of variation of support damping diminishes. 

5. The off-tuned support (K i M) can be designed to produce a 
consideraDle improvement in system response in comparison to the rotor on 
rigid supports. If insufficient damping is incorporated in The supper 
then the resulting rotor steady-state amplitude may be larger than the 
original rotor response for support rJ Mffness values K > I. 

6. If there is excessive suppor + damping (C > 20) with a low mass 
ratio support (M = 0.1), then the forces transmitted through the support 
may exceed the unbalance forces (TRDS > 1.0). 

7. Although the steady-state analysis shows that the rotor amplitude 
will be small for an underdamped (C < 0.50) low mass ratio support system, 
the orbital analysis shows that a large initial transient motion can be 
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generated due to the suddenly applied unbalance force and that this motion 
is not readily attenuated. 


8. The optimum damping based on minimization of the rotor steady- 
state amplitude for both tuned and off-tuned conditions produces a satis- 
factory transient response from the standpoint of rapid reduction of the 
initial transient motion, improved system stability and reduction of the 
forces transmitted. 
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APPENDIX A 
DISCUSSION 

A high speed rotor shaft may be considered as a continuous elastic 
member with variable mass and inertia properties along its length. The 
rotor shaft usually has attached to it such components as turbine or 
compressor blades, impeller disks, or spacer assemblies or seals. If 
the axial dimensions of each rotor component is small in comparison to the 
overal I length of the rotor, then each compgnen+ may be treated as a con- 
centrated mass with a polar moment of inertia equivalent to that of the 
original component. If the mass of the components are large in comparison 
to the shaft mass connecting the components, then the shaft weight can 
be neglected or considered to be located at the mass stations. If fne polar 
moment of inertia of each section is ignored, then the stations may be 
considered as concentrated masses, rather than distributed in the plane 
of the rotor element. However, if the sections whirl in a plane, perpen- 
dicular to the spin axis then the gyroscopic moments do not act on the 
system and hence the equations reduce to the same as if point masses were 
assumed. 

The position vector of the nth mass center is given by 

P = T + 7. + 7 + "e 

b j s u 


where 

6^ = vectoral bearing support deflection 

= vectoral journal deflection 

6 = vectoral shaft deflection 

s 

e* u = displacement of mass center from the shaft centerline 
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The system being analyzed has been reduced to a single mass rotor 
mounted in idealized linear bearings and the bearings are in turn mounted 
on damped, elastic supports. By considering only small deflections, the 
spring rate of the flexible, massless rotor shaft may be considered to be 
I inear. 

The rotor disk (see Fig. I) is considered to whirl in a plane and 
hence no gyroscopic moments are acting on the system. The orthogonal support 
and bearing spring rates are assumed symmetric and no cross coupling terms 
are considered to be acting at the support housings. The aforementioned 
assumptions al lov* me equations of motion of the system to be written as 
total differential equations- 

DERIVATION OF EQUATIONS OF MOTION 


A. I Kinematics 

The position vectors to the mass stations are given by 


m l : bearing Mi/0 

housing mass P = + Y i^y 

m 2 : rotor mass 


(A. I ) 


M 2 /0 

= (X 2 + e cos0)Ti + (Yo + e s I n© )7? 
z u x ^ u y 


The velocities of the mass stations are given by 

Mj/0 . . 

V = XiirT + YirT 

‘x 1 y 

_M ? /0 . . 

"V = (X 2 - e 0sin0)‘n + (Y 2 + e 0cos0)"n" 
*■ u x * u y 


(A. 2) 


(A. 3) 
(A. 4 ) 
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A. 2 Kinetic Energy 

The kinetic energy of the system is given by 

<Vi • > + r i, . 

i=i j 

neglecting the gyroscopic couples acting on the disc, the system kinetic 
energy reduces to 

• • 

T * 4 * .^{Xi 2 + Y! 2 ) + 4 M 2 [(X 2 - e 0sin0) 2 + (y 2 + e ,0cos0) 2 .l 

Z 4 U J-l 

+ i $ 0 2 (A. 6) 

2 zz 

A. 3 Potential Energy 

The potential energy of the system is composed of the sums of the 
potential energy of the flexible shaft, the potential energy of the bearings, 
and the energy of support structure as fol lows 

V = i [K (X 2 + Y 2 ) + K (X. 2 t r . 2 ) h MX! 2 1 ' 1 2 )J (A. 7) 

/ S S S b J J 



3 

7 <J>. .w.6 j. (A. 5) 

=1 *J 1 J 


where 


x s = x 2 - Xj - X. 

A. 4 Dissipative Energy 

The system dissipative energy consists of the damping functions provided 
by the bearing support system, the bearings, and the external and internal 
rotor damping and the aerodynamic rotor cross coupling as described by 
Alford (4). 

D = ~ { C j (X ] 2 + Yj 2 ) + C b (Xj 2 + Yj 2 ) + C.(X 2 2 + Y 2 2 ) 

+ C.[x 2 + Y 2 + 2 oj (Y X - X Y )] + Q(Y 2 X 2 - X 2 Y 2 )} (A. 8 ) 

l u S S SS SS ^ ^ ^ 4 
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The internal damping function is dependent upon the rotor precession rate 
and can cause self excited whirl instability when the rotor is operated 
above the c-ilical speed (26). Alford has demonstrated that the aerodynamic 
cross coupling stiffness term can also cause -otor instability when the 
rotor speed is supercritical . When the system dissipation function is com- 
prised of only the first three terms, the system is inherently stable. 


A. 5 Lagranges Equations . 

The governing equations of motion are obtained from Lagranges Equations 
which state: 


d+ 3A 
3q r 


3L + 3D_ _ p 
3q r 3q‘ r 


(A. 9 ) 


where 


L = T - V 

The total number of equations of motion obtained will be equal to the 
number of degrees of freedom of the system which is seven and are given as 
fol lows. 


Rotor 


< 2 : M 2 X 2 + C s X 2 + C.(X 2 - Xj - X.) + K s (X 2 - X L - X.) 

+ QY 2 + wC.(Y 2 - Yi - Y J = M 2 e u (w 2 cos(ot) + osin(ut)) 


(A. 10) 


Y 2 : M 2 Y 2 +CY 2 +C.(Y 2 -V 1 • Y.) x K (Y 2 - Yj - Y.) - QX 2 

Si J *• J 


(A. I I ) 


- u)C.(X 2 - Xj - . > = Mge u (u) 2 sin(<i>t) - asin(vf): 
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Bear! ngs 


X. : (C. + CJX. - C.(X 2 - Xi) + (K + K )X. 

J i b j i z 1 b s j 

- K s (X 2 - X!) - C.w(Y 2 - Yj - Yj) = 0 

Y. : (C. + C.)Y. - C.(Y 2 - Y , ) + (K. + K )Y. 

J i b j i z b s j 

- K s (Y 2 - Yx ) + C.o)(X 2 - Xj - X.) = 0 

Support 

Xj: MiXj + (Ci + C. )Xj - Cj(X 2 - X. ) + (Kj + K^Xj 

- K (X 2 - X.) - C.w(Y 2 - Y, - Y.) = A 

s z J i z 1 J 

Y j : MjYj + (Cj + C. )Yj - C. (Y 2 - Yj ) + (Kj + K g )'» 1 

- K CY 2 - Y.) + C.ai(X 2 - Xi - X.) ■= 0 

s ^ J \ 1 1 j 

Angular Acceleration 

0: ^ zz + M 2 e ^ ) 6 + M 2 e|[Y 2 cos C - X 2 sin0 

• • * * 

- 0iY 2 s’n0 + X 2 cos0)] = T (6 ) 


i A . 12) 


(A. 13' 


(A. 14) 


(A. 15) 


(A. 16) 


Where 

0 = u>, 0 = u 

The equations A. 10 to A. 15 may be vectorially combined by representing 
the displacements in complex notation as fellows 
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Z 2 » x 2 + iY 2 

Z. = X; + iY. 

Zi = Xj + iY 4 

; t ^ 4. r (7~ - Zt -Z-) + K.(Z 2 - Zj - Z.) 

4 2 : i*<2 £ -2 * J ■* 

- iQZ 2 - i«Cj(Z 2 Zj - Zj) = M 2 e u (« 2 - ia)e'“ + 

Z j : (c b + c.)i. - c.a 2 - ii) + ^ + K s >z j 

- K s (Z 2 - Zi ) + SC.w(Z 2 - Zi - Zj) ■ o 

Zl : + (C! + C.JZ! - C.(Z 2 - Zj) + («<i + * s >Zl 

- K (Z 2 - Z.) + iC.u(Z 2 - Zi - Zj) » 0 

s *■ J 1 J 


(a.:7) 


CA.I3) 


(A. 19) 


(A.20) 


w . 


35 












■ -Yi.Tttjn-fcriftwil I, i 




I 


*&-•? . *t><\,'-> '■* yvs<5r r *‘ : f p| * • 


I ___ 

4 .„ — -« «»r*4s. WUWJ'JMMJBSI ! 


APPENDIX B 

TRANSIENT COMPUTER PROGRAM FOR THE THREE MASS SYSTEM 

The rotor system used in the transient analysis is of considerably 
greater -~omplexity than the equations used in the steady-state analysis, 
in the more general transient analysis, internal and aerodynamic cross 
coupling may be incorporated in the rotor. The bearing may have 8 stiff- 
ness and damping coefficients while the support system may have 4 stiffness 
and damping coefficients. Because of the more generalized treatment, 
studies may be conducted on the stability of the rotor due to hydrodynamic 
bearing forces, internal friction, or aerodynamic cross coupling. The 
support system can be investigated to show the influence of the support- 
damping in promoting stability. 

Equations of Motion 

X 2 : M 2 X 2 + C g X 2 + C. (X 2 - Xj - Xj) 

+ K s (X 2 - Xj - X.) + QY 2 + uC.(Y 2 - Yj - YJ 

(e. i) 


+ K s (Y 2 - Yi - Yj> - QX 2 - wC.(X 2 - X l - Xj) 

= M 2 e u &» 2 sina>t (B.2) 

X.: 2M.(X, + X.) + (C. + 2C )X. - C. (X 2 - X 2 ) 

J J J « xx j i £ 1 

+ (2K + K )X. - K (X 2 - Xi) - C.u(Y 2 - Y, - Y.) 

XX S J S z 1 I J 

+ 2(K Y. + C Y ) = 0 (B.3) 

xy j xy 


= M 2 e u w 2 cosut 


Y 2 : M 2 Y 2 + C s Y 2 + C.(Y 2 - Y x 
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(B.4) 


Y.: 

J 


+ Yj) + (C. + 2C yy >Yj - C.(Y 2 - Y x ) 


(2K 


yy 


+ K )Y - K (Y 2 - Yj) + C.u(X 2 - Xj - X.) 

5 J S I J 


+ 2(K X. + C X.) = 0 
yx j yx j 


X x : 2MXX! + 2M.{Xi + X.) + (2C + C. )X x 

J J lx i 

- C.(X 2 - X.) + (2K + K )X» - K (X 2 - X.) 

i 2 j lx s 1 s z j 

- C.u(Y 2 - Yi - Yj) ■ 0 

Y x : 2M 1 Y 1 + 2M. (Yi + Y.) + (2C, + C. )Y, 

1 11 J J iy i 4 

- C. (Y 2 - Y.) + (2K + K )Yi - K (Y 2 - Y.) 

i 2 j iy s 1 s 2 j 


+ C.w(X 2 - Xi - X.) = 0 

I J 


The support equations can be reduced to the following 
Eqs. B.3 and B.4. 


MjXj + C K) + K Xx - C X. - C Y. 
11 lx 1 ix 1 xx j xy j 


- K X. - K Y. = 0 
xx j xy j 


M 1 Y 1 + C ,Yi + K Yi - C X. - C Y. 
11 iy 1 iy 1 yx j yy j 


- K X. - K Y. = 0 

yx j yy j 
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(B.5) 


(B.6) 

form by using 


(B. 7) 


(B.8) 
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11*55 AM JULY 09* 197} ***************** 

BEGIN 

COMMENT 

THIS program CALCULATES the TRANSIENT RESPONSE or A SINGLE-MASS 
FLEXIBLE ROTOR ON FLEXIBLE BEARING SUPPORTS (LINEAR) * THE RESULTS 
ARE PLOTTED AUTOMATICALLY (SEE INPUT CAROS BELOH) AMO THE 
INITIAL ANO FINAL VALUES ARE PRINTED OUT ON THE LINE PRINTER, 

THEY APPEAR In THE CORRECT ORDER Tu BE TYPED ON DATA CAROS 
TO CONTINUE THE MOTION IF DESIRED (THAT IS CARD 5(A)), 

END OF COMMENT! 


the input to the program is as follows i 

(ALL OATA IS IN FREE FIELD ) 

comment 

card i reao;cr«/»tmax*h*n*casno) 

T MAX - NO, OF RADIANS SOL, IS TO BE SOLVEO 
H • STEP SIZE USEO IN INTEGRATION PROCEDURE ( TMAX/HS700) 

N * NO* OF EOS, TO BE INTFGRATEO 

( • 12 IF SUPPORT • • B IF NO SUPPORT) 

CASNO - IDENTIFICATION NUMBER (XXXX,XXXX) 

t(MQ,)(OAY),(YEAR)(CASE NO,)] 

CARD 2 READ(CR*/»HPN»M*K*OS*Cl*QAC»EU.Nj»Nl) 

RPM • ROTOR SPEED * REV/MIN 
H - ROTOR HEIGHT * LB* 

K * ROTOR SHAFT STIFFNESS * LB/IN 

OS - ABSOLUTE SHAFT DAMPING * LB-SEC/IN 

C! • INTERNAL FRICTION DAMPING * LB-SEC/IN 

0 AC * CROSS-COUPLING* LB/tN 

EU - UNBALANCE ECCENTRICITY OF M IN MILS 

mmmm 

HJ - JOURNAL HEIGHT AT EACH ENO * LB. 

Hi - SUPPORT HEIGHT AT EACH BEARING * LB, 

C»RO 3 READ(CR#/*KXX*KYY*CXX»CYY»KXY»KYX*CXY#CYX> 

BEARING STIFFNESS AHO DAMPING of EACH BEARING 

CARO 4 IF N»B THEN REAO(CR*/.KlX*KlY#CtX*ClY) 

SUPPORT STIFFNESS AND DAMPING OF EACH SUPPORT 

C*RO 5 (4) READ(CR»/.FOR 1*0 STEP 1 UNTIL N DO(YCO»in) 

0 IS FOR INITIAL TINE* RADIANS 

1 - ABSOLUTE ROTOR OISP. * X-OIR, 

2 - ABSOLUTE ROTOR VELOCITY* X-DIK, 

3 - ABS, ROTOR OISP, * T-OIR, 

4 - ABS, ROTOR VEL. * Y-DIR, 

5*6*7»6 - SAME ORDER AS ABOVE FOR JOURNAL RELATIVE NOTION 
9*10*11*12 - SAME AS ABOVE FOR SUPPORT MOTION (IF REOUIREO) 

CARD 6 (5)REA0(CR»/CS» ISCALE*XHIN»0X*XMIN2#0X2*XMIN3»0X3) 

CS • PLOTTER CONTROL 0 - NO PLOT 1 • PLOT 

ISCALE • SCALE CONTROL 0 - PROG, SCALE 1 - USE FOLLONING INFO, 
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XMIN - NO. TO APPEAR AT ORIGIN ROTOR OISLPACEMENT 

OX • SCALE INCREMENT PER INCH ROTOR DISPLACEMENT 

XMIN2#0X2 • SAME AS ABOVE BUT FOR JOURNAL PLOT 
XM |N3#DX J • SAME AS ABOVE BUT FOR SUPPORT 

CARD 7 (6) IF CSAO THEN REAd<CR#S*RPCS) 

RPCS • NO. Or TIMES INTEG. FOR TMAX RAOIANS IS To BE SOLVEO 

example DATA INPUT I 
12,56,0.05*6*606.7101* 
lOOOO* 675# 260000# 0*0# 0*0, 5* 312*50# 

351000*606000*739,665* <>,0,0,0, 

0#»0 ,5*0,0# -0 ,5#*0 .05# 0*0* *0,05* 

1 * 0 # 0 , 0 , 0 * 0 # 0 , 0 , 

5, 


NOTE ON PLOTTING I 

REQUEST ONE (1) BLOCK PER SET OF DATA 
END OF COMMENT I 

X •-•REPEAT THE SERIES FOR EACH CASE — 

X 

INTEGER CS I 
INTEGER I*J»NA,N) 
boolean I SCALE I 
boolean stable I 

INTEGER RPC#RPCS| 

INTEGER VV 1 

ALPHA ARRAY ALPl# *LP2tOI61, ALP3 # AlPA# ALPS, AlPO# ALP7, AlP 6, 
ALP9»ALPlO,ALPll,ALPl2*ALPl3,ALPU,ALPl5#ALPl6»ALPl7tO«3II 
ALPHA ARRAY ALPHA!, ALPHA2(0«2Jj 
ALPHA ARRAY ALP16 , ALP1 9(0 1 A )| 

ALPHA ARRAY ALP20(0t61* ALP2ltO>6], ALP22(0»3U 
ALPHA ARRAY ALP23*ALP2AC0l«I| 

ALPHA ARRAY ALP25*ALP26#ALP27*ALP28(0I71) 

ALPHA ARRAY ALP29* ALP30 (0 1 31 1 
REAL CASNOI 

REAL K2X*K2Y,C2X,C2Y, MCX, MCY* ACX# ACY ) 

REAL K11,K12,K13,K14,FMBH,FMBI*FMSH,FMSI,TIMHB,TIMHSI 
REAL TMAXH | 

REAL FMSHH,FMBHH,TIMH8H,TTMHSHJ 

REAL 0AC,yMN2,DX2*XMAX2,XM!N3#DX3,XHAX3, YMIN*,DY2, YMAX2* 
YMIN3*DV3*YMAX3I 

REAL ZZ#ZZ2*Z1#Z2»23,Z4#Z5#Z6*Z7,Z6,Z9»Z10 ,Y0»Y1, Y2# Y3#Y«,Y5. 

Y6,Y7#Y6#Y9#Y|0J 

REAL Y11*Y12,DS*C1*MJ,KXX*KYY*CXX#CYY*K!X*K1Y»C1X*C1YI 
REAL MJ,Z22,Z32#ZU,Z12*Z13*Zl«*Z15*Z16*Z17*Zld ,EU ) 

REAL XA,HXA*XI,HXI*YA#HYA#YI,HY1 I 
REAL TMAX , H,CC, R, OMEGA, 0MEGA2 I 
REAL NCS I 

REAL DUM1,OUM2,DUN3,DIH4,OUM22,OUM42,RE1»IM1,M1.N1,C!,K1*KB,CB I 
REAL K2,C2 I 

REAL RPM,N,K,CD#CQ,DC,RAD,M,00C,NC1,NCRATI0«DEN1,AA,0EN2,ACR, 
XMAX, YMAX \ 

REAL XMIN# YMIN# OX, OY I 
REAL KXY,KYX#CXY,CYX} 

REAL Z19,Z20,Z21,Z23) 
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LABEL ALLDONE # ACARD I 

REAL ARRAY A,B,C tO I A1#«#KK# Y {01 4*01 121 #F (0« 12)# 
ARRAY HXR#BYR[0l700)! 

LABEL REPEAT! 

FORMAT OUTDATA (X2*EI3.6#X2»F6.2»4(X2#E13»6))! 


PROCEDURE MAXMIN <X,Y#N*M)! VALUE N»N! 

REAL HJ INTEGER N! ARRAY X»YlOI I 

BEGIN 

INTEGER I I 

XI ♦ XA ♦ Xti) | VI ♦ VA ♦ Ytl) | 

FOR I ♦ 2 STEP I UNTIL N DO 
BEGIN 

IF XtlJ > XA THEN 

begin 

XA » X[II! MX A * I | 

END 

else begin 

IF XtU * XI THEN BEGIN XI ♦ XIII ! Hxl * I I END 
IF YU) > YA ThEN 
BEGIN 

YA ♦ YIDI HYA ♦ I ) 

END 

ELSE BEGIN IF Ytl J < VI THEN BEGIN Yl ♦ YCI) ! HYI 
END ! 

HYA ♦ (HYA-I) x H / 6.28 | 

HYI «■ (HYI-1) x H / 6.28 I 

HX I ♦ (HXI-1 ) x H / 6.28 I 

HXA ♦ ( HXA» 1 ) x H / 6.28 ! 

END DF MAXMIN ! 

prccldurf samescale(x#y#n»xnin>xmax*dx»ymin*ymax#dy)! 

ARRAY X.YIO) I INTEGER N I 
REAL XMIN»XHAX»DX»VMIN#YMAX*DY! 

begin 

REAL ANS, ANSX ) 

INTEGER l| 

ANS ♦ ABS (YU))! AN$X * ADS (XU))! 

FOR I «• 2 STEP 1 UNTIL N DO 
BEGIN 

IF ABS (YII)) > ANS THEN ANS ♦ ABS (YCI))) 

IF APS (XII ) ) > ANSX THEN ANSX ♦ ABS (XtJ)>! 

END! 

IF ANSX > ANS THEN ANS ♦ ANSX! 

IF ANS < 1 THEN 
8EGIN 

XMIN ♦ *11 DX ♦ 0*3333333333! 

END 

ELSE 

IF ANS < 3 THEN 
BEGIN 

XMIN «• -3 J DX ♦ I ! 

END 

ELSE 

IF ANS < 6 THEN 


12*01700] I 


END | 


* I ! END! END! 


VALUE N| 


BEGIN 

XNIN ♦ *61 
DX * 21 
END 

ELSE IF ANS <12 THEN 
BEGIN 

XMIN **121 
OX *41 
CNO 

ELSE IF ANS < 18 THEN 
BEGIN 

XMIN ♦ M8 I OX ♦ 6 | 

END 

ELSE 

BEGIN 

XMJN ♦ - 24 J 
OX *8j 
ENQI 

OY ♦ OX; YMIN * XMIN| 

XMAX ♦ 3*0X1 YMAX ♦ 3*0Y| 

END of sanescale I 

PROCEDURE PLOTCHEK(X#Y#N#XMIN#XMAX#YMIN#YMAX)l VALUE Nl 
ARRAY X,Y tO J ! INTEGER Nl 
REAL XMIN# XMAX# YMIN# YHAXI 
BEGIN 

INTEGER II 

FOR I * l STEP 1 UNTIL N 00 
BEGIN 

IF X til > XMAX THEN X (II ♦ XMAX 
ELSE 

IF XIII < XMIN THEN XT I J ♦ XMIN! 

IF Y II) > YMAX THEN Y(I] ♦ YMAX 
ELSE 

IF VIII < YMIN THEN Y 1 1 J * YM IN | 

END I 

END OF PLOTCHERI 


PROCEDURE ORBITTOP (1)1 
VALUE II INTEGER I I 
BEGIN 

IF I • 1 THEN 

SYNBOLC 1 ,50# 9 ,00# #21# ALP1 #0# 24 ) 

ELSE IF I » 2 THEN 

SYMBOL(2,00#9.00##21#ALP2#0#18) 

ELSE 

SYMB0L(2.00#9.00# # 21#ALP3#0#16)I 
SYMBOL(5,00#9.50#C,10#ALP30#0#3)I 
NUMBER (5. 20# 9, 50# 0. 10#CASN0#0»4)I 
SYNBOLC 1,55#8.T5#» 14# ALP25#0#2G)I 
NUMBER(2,03#B.75#.U#AYC0#1)/6.28»0#2)I 

K ♦K/1000J KXX«-KXX# 10001 KYY*KW/ 1000! K 1X*K IX/ 1000IK 1 Y*K IY/ 1000| 
SYMB0L(0,75#6,50#,14.ALP4 #0# 18 )| NUMBER ( 1 <35 #6 *50# * 14 #RPM #0,0)1 
SYMB0L(3, 75,8.50# ( 14#ALP6 #0# 19) | NUMBER (4 ,23# 8 .50# , 14# EU #0,3)1 
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SYMBOL 0,75*6,23.* ♦4*^>U , S#Q»UHMUMBCBC1 ,2M ,25*,M*N 


Ii2$>,U>NJ 
0,00* ,14#K 
0,00*,14#KXX 
7 *75# ,14#DS 
7,75».14#KYY 
7,50#,U#CI 
7*50# ,1*»CXX 
7,25#,14#0AC 
7 *25# ,14#CYY 
7 t 00*,14*MCX 
7*00# ,|4»ACX 
0*75#,|4#NCY 
6,75#,1A*ACY 

6 *50# * 1* #Ml 
6*27# «|4#K1X 
4,27#.14#C1X 
4 *05# , i* #KI Y 
4*05# ,14#C1Y 


*0*2)1 
# 0 # 2 )) 
#0# 3) I 
#0*3)1 
# 0 * 2)1 
#0*3)1 
# 0 * 2)1 
# 0 * 2)1 
# 0 * 2 )) 
* 0 * 2)1 
» 0 * 2 )f 
* 0 * 2)1 
* 0 « 2 )| 
* 0 * 2 )) 


SYMBOL! 3*75*6 ,25#, 14* AI.P6 #0*17) JNUH0ER<4,23 
SYMBOL (0. 75*6,00#, 14»ALP7 #0#21) )NUMBER!1,23 
SYMBOL! 3, 75*6 ,00*, L4#ALPl0*0*22)fNUM6ER(4,35 
0YPBOL(0,75#7,75#,14#ALP9 • 0# 23) I NUMBER (1 ,23 
SYMB0U3,75#7,75#,14»ALP12#0»22)|NUMBER(4,35 
SYMBOKO, 75*7*50# , 1 4#Al.Pll #0# 23))NUMBER( | , 23 
SYMB0U3,75#7,30#,14#ALP14#0,24))NUMBER(4.35 
SYMBOL!©, 75*7, 25*. 14#ALP13#0#16 )jnuMBER! 1, 12 
SYMB0U3, 75*7,25#, 14* ALPt4#0#24))NUMBER(4. 35 
SYMBOKO ,75*7. 00 »,14»A|.P15*0*t9)|NUMBER( 1.47 
SYMBOL (3 *75*7, 00#, 14*ALP18#0#4 ))NL.4BER!4,23 
SYMBOKO, 75*6, 75* ,14* ALP17*0> 19 )| nUMBERC 1, 47 
SYMBOL (3, 75* 6, 75# ,14* ALP20#0#6 )|NUMBE*!4,23 
IF N > 6 THEM BEGIN 
S YMBOLIO, 75*4, 50*. 14*ALP19*0*16 )|NUMBER( 1,23*4, 50#, 14 *N1 #0*2) f 

SYMBOL (0 ,75*4,27*, 14# ALP21#0, 21 )|NUMBERC 1,23*6, 27# ,14#K1X #0*3)1 

SYMBOL! 3,75*4, 27#, 14* ALF22*0*24 ))NUMBER(4, 35*4, 27*, 14#C1X #0*2)1 

SYMBOLS. 75*6, 05#, 14 *ALP23»0*21)fNUMBEHU ,23*4 ,05* , 14,K1 Y *0*3)1 

SYMBOL(3,75*4,05#.14*ALP24*0*24))NUMBER(4.35*4,OS*,14*C1Y *0*2)) 

END) 

K ♦XKlOOOl KXX*KXX*1000I KYY,KYYM000) K1 X*KIXk 100O)K1Y*K1Y*1OOOJ 
END Of ORBITTOP I 

PROCEDURE AGRID ! XMlN#DX#YNlN»OY)) 

REAL XMIN#OX*YHIN«DYj 
BEGIN 

AX IS ( 0,0 #0# ALPHA 1* 14#4#90*YHlN#0Y)) 
AXIS(0#0#ALPHA2*-14*6»0»XMIH,OX)| 

AXI S(6* 0* ALPHA 1# 0*6# 90* YMIN,DY)J 
AXlS(4#6*ALPHA2*0*6*ie0#XMIN*0X)| 

END# 

PROCEDURE FUNCTION <J#Y)I 

VALUE Jl INTEGER J) ARRAY Yt0#01) 

BEGIN 

YO*YCJ#01) 

Yl*VU#ll) 

V2*Y!J#21) 

Y3+YIJ# 31 1 
Y4*Yt J#41) 

Y5*Y!J#51) 

Y6*YIJ#6)I 

Y7»YCJ#71I 

YB»YIJ*B)) 

IF N * 6 THEN BEGIN 
Y9*Yt J#9) I 
V10 ♦ YCJ#10J) 

Yll «• Y I J# 1 1 ) I 
V12 ♦ Y( J# 12) I 
END) 

f 101 ♦ l) 
fill ♦ Y2! 

F12J ♦ EUkC0S(Y0)-Z1*Y2»72*(Y2-Y10-Y4) 

"Z3*C Y1»Y9»Y5 )*24* Y3»72*< Y3*Yl l*Y7 )) 

F 1 31 * Y4) 

F14 1 ♦ EUkSIN(Y0)»Z1kY4*Z2k<Y4»Y12»Y8) 


42 


i 



I 


•Z3x(Y3*Ylt*Y7)4Z«*Yl4Z2*<Yt*Y9 -Y5)| 

IF N > 6 THEN begin 
FC9) ♦ Y10| 

FUOl 4 •ZiOxcZllxY94Zl2*VlO*Z13xY5-ZtA*Y6 -Z l9xY7«Z?0xYe) J 

run * Yi 2 i 

Ft 12 1 ♦ *Z10*<Zt5KYU*Z«*Y12-Zl7XYWl»>«Ya*Z21*V5-Z23*Y6}! 
END ELSE Ft9J*niOJ*Ff tlJ*Ft|21*0.0| 

Ft53 ♦ Y6I 

F161 ♦ -Ftl01-Z5x(Z6xY6-Z22x(Y2»Y10> 4Z7xY5 
•Z32xCYl«Y9)«Z22x<Y3*Yn-Yr>*Z19xYZ*Z20xY§)l 
Ft71 ♦ Y0I 

FtGJ ♦ -Ftl2J -Z5*(Z«xYa-722x(Y4»Y12) 4 Z9xY7 
•Z32XCY3-Y11 ) 4 Z22x( Y1*Y9*Y5 )*Z21*Y54Z23* Y6 ) J 

END OF FUNCTION I 
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PROCEDURE RK6 <K/ Y) j INTEGER K l 

ARRAY Y EO# 0 3 I 

BEGIN 

REAL P I INTEGER I*J I 

FOR J ♦ 1 STEP 1 UNTIL A 00 

BEGIN 

FUNCTION CJ-l/Y) I 

FOR I ♦ 0 STEP 1 UNTIL N 00 

KKCJ/I) 4 fill | 

FOR I 4 0 STEP l UNTIL N 00 
BEGIN 

P 4 At J J x (KKtJ,U - BtJJ x QtJ-1,1)) ) 

YtJ/II 4 YrJ-l.IJ 4 H x P | 

Otj/I) ♦ Of J#l/I) ♦ 3 x P - Ctjl * KKt J.I1 I 

ENDI 

ENDl 

FOR 140 STEP 1 UNTIL N 00 

begin 

VtO>IJ4YtA#IJI 

oto#i)4at«#n i 

AY(I#KJ 4 YIA#UJ 

ENOJ 

FNBI *(Zl3xYfO/5J 4 Zl«xYfO/6J 4 Z|9xYtO/7J 4 Z20xY tO/8J )*2 

4 (Z21xY 10#5J 4 723XYI 0/6 J 4 Zl7xYtO/7J 4 ZlBxYtO/6 J >*2 I 
FH SI 4 (Zll*YtO/9l4Zi2xYrC/lO])424<Zl5xYtO/llJ4Zl6xYtO/i2J>*2 > 
IF FMBI>FMBH THEN BEGIN 
FMBH 4 FNBI I TIMH8 4 Y(0«01J ENOj 
IF FMSI > FMSH THEN BEGIN 
FMSH 4 FMSI I TIMMS 4 YtO»Ol| ENOJ 
ENO OF RKGl 

PROCEDURE TIMESTEP ( TMAX/ H/ N / AY/ NA/YIJ 

VALUE TMAX/H/NJ 

REAL TMAX/H! INTEGER N/ NA J 

REAL ARRAY AYIO/JJ ) 

ARRAY YtO/Ol I 
BEGIN 

INTEGER I #J/M| LABEL REPEAT J 
FOR 140 STEP l UNTIL N 00 

BEGIN 


43 



i a i « 5 1 t , - . igr». - ■ if ^= * 


END! 


6(0#!1«-0) 

AY[i # n*rto#ni 


i 


FHBH ♦CZ13*YC0»51 ♦ Z14*YfO#6) ♦ Z19KYtO#7] ♦ Z2d*YlO#8] >*2 

♦ (Z21*Y[0*51 ♦ Z23*Y f 0* 6) ♦ Z17*Y[0,7J ♦ ZIB*Y10#8) )*;> J 

FMSH ♦ (Zil*Yt0#9)*Zl2KYt0#!0n*2*(Zi5>«Yr 0# 1 1 )*Zl 6*YI 0* 12 1 )»2 I 
TIMHB ♦ Yl0#01j TJMHS ♦ Y(0#0)| 

KMJ 

REPEAT* 

K+K+1J 
RKG (K#Y )l 

IF YtO#OKTMA?: THEN GO TO REPEAT I 
N A*K j 

ENO OF TIMESTEPJ 


Atl] «• CI1J * CI4 ] ♦ 0,5 | 
At 21 * CC 23 «• 1-S0RT{0,5) | 
At 31 ♦ Ct 33 * i*SBRT (0,5) | 
At 4 3 ♦ 1/6 I 
etl] ♦ BE* 3 ♦ 2 ; 


Bt 2 1 «■ BE 3 3 ♦ II 

HRITr(LP#< / /. X1Q*«****** SMFROLS ******".//» 

" THIS PROGRAM CALCULATES THE TRANSIENT RESPONSE OF A SINGLE-MASS"#/# 
"FLEXIBLE ROTOR ON FLEXIBLE BEARING SUPPORTS (LINEAR) , THE RESULTS"# 

/#* are plotted automatically (see input caros below) and the •*/# 

"INITIAL AND FINAL VALUES ARE PRINTED OUT ON THE LINE PRINTER. "#/# 

"THEY APPEAR IN THE CORRECT ORDER TO BE TYPED ON DATA CARDS "*/# 

"TO CONTINUE THE MOTION IF DESIRED (THAT IS CARD 5<A>>, "#/# 

///»" THE INPUT TO THE PROGRAM IS AS FOLLOWS I "*/< 

" (ALL DATA IS IN FREE FIELD ) "»//# 

" CARD 1 RE A0( CR#/, TMAX,H,N*CASNO) "»/» 

" T MAX - NO. OF RADIANS SOL. IS TO BE SOLVED *»/* 

" H - STEP SIZE USED IN INTEGRATION PROCEDURE ( TNAX/HS7 00) "»/» 

" N - NO, OF EOS, TO BE INTEGRATED "»/* 

" ( « 12 IF SUPPORT * » 8 IF NO SUPPORT) "*/# 

" CASNO - IDENTIFICATION NUMBER tXXXX.XXXX) "*/* 

" l(MO,)(OAY),(YEAK)(CASE NO.)) **/* 

//," CARD 2 READ(CR>/,RPM#W»K»DS,CI*GAC#EU»NJ*N1) "#/» 

" RPM - ROTOR SPEED # REV/MIN "»/# 


ft 


ft 


//." 


//." 


//." 

m 


M • ROTOR WEIGHT , LB. "./# 

K - ROTOR SHAFT STIFFNESS » LB/IN 

DS - ABSOLUTE SHAFT DAMPING * LB-SEC/IN 

Cl - INTERNAL FRICTION DAMPING * LB-SEC/IN 

OAC • CROSS-COUPLING, LB/IN 

EU • UNBALANCE ECCENTRICITY OF N IN MILS 


*,/# 


"#/, 

*»/» 


WJ • JOURNAL WEIGHT AT EACH END , LB* "»/# 

W1 • SUPPORT WEIGHT AT EACH BEARING # LB, "*/# 

CARD 3 READ(CR,/.KXX#KYY*CXX*CYY»KXY»KYX»CXV#CYX) "»/# 

BEARING STIFFNESS AND DAMPING OF EACH BEARING "»/# 

CARD A IF N>6 THEN RE AD(CR»/* MX# M Y» C1X* Cl Y ) "*/.* 

SUPPORT STIFFNESS AND DAMPING OF EACH SUPPORT 


CARD 5 (4) REAO(CR»/»FOR 1*0 STEP 1 UNTIL N DOE YEO* 1333 
0 IS FOR INITIAL TIME* RADIANS "»/# 


"#/» 
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J. jLX 1 1 i ii*i 1, . U a JL . i JkA 4*** 4- ..i.'i M. A .Jlu 


i 


* 1 • ABSOLUTE rotor oisp, * X-OIR, "»/, 

" 2 • ABSOLUTE ROTOR VELOCITY# X"DIK. ",/# 

" 3 • ABS, ROTOR OISP, # Y-DIR, ",/, 

" A • ABS, ROTOR VEL# # Y"DIR, "»/, 

* 5*6,7, 6 • SAME OROER AS ABOVE FOR JOURNAL RELATIVE MOTION ",/, 


" 9# 10# 1 1# 12 - SAME AS ABOVE FOR SUPPORT MOTION (IF REQUIRED) " 

//#" CARO 6 <5)READ(CR,/CS#ISCALE,xMIN,DX»XMIN2,Dx2*XMlN3,0x3) 

* C5 - PLOTTER CONTROL 0 - NO PLOT 1 • PLOT "»/# 

" ISC ALL • SCALE CONTROL 0 - PROG. SCALE 1 • USE FOLLOWING INFO," 

#/.• XMI N - NO. TO APPEAR AT ORIGIN ROTOR OISLPACEMENT •#/# 

" OX • SCALE INCREMENT PER INCH ROTOR DISPLACEMENT ",/, 

" XMIn?#DX2 • SAME AS ABOVE BUT FOR JOURNAL PLOT "#/# 

* XMIN3#0X3 • SAME AS ABOVE BUT FOR SUPPORT "#/# 

//," CARO 7 (6) IF CSpO THEN REAOCCR,/,RPCS> ",/, 


" RPCS - NO. OF TIMES INTEG, FOR TMAX RAOIANS IS TO BE SOLVED "#/» 
//#" EXAMPLE OATA INPUTI ",/# 

"12.56#0.05#8#608.7t01# "#/# 

" 10000# 675# 280000# 0»Q#0# 0,5# 312* 30# "#/# 

"35 1000 #60 6000, 7 39, 865# ",/r 

" 0 # 0 , 0 , 0 , 0 , 0 , 0 , 0 # 0 , ",/# 

" 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , "#/, 

"5, "#/, 

//#" NOTE ON PLOTTING! "#/, 

" REQUEST ONE (I) BLOCK PER SET OF DATA ",/, 

" —REPEAT THE SERIES FOR EACH CASE — ">)l 

WRITE ( LPI PAGE I ) I 

NRITF(LP»<" PLOTTER OUTPUT INFORMATION AND SUGGESTIONS",//. 
"MAKE TMAX A MULTIPLE OF 6,?® (BUT LESS THAN ’5 IF H«Q.05> ",/, 

"A REASONABLE VALUE OF H IS 0,05 WHICH GIVES 125 STEPS PER "#/» 
"CYCLE OF RUNNING SPEED FOR THE INTEGRATION,",//# 

"TWO (?) MINUTES PROCESSOR TIME IS REQUIRED FOR 10 CYCLES OF ",/. 
"SOLUTION FOR N « 6, ",/, 

"THREE (3) MINUTES PROCESSOR TIME Is REQUIRED FOR JO CYCLES OF »,/, 
"SOLUTION FOR N ■ 12, "»/# 

/# "A SMALL CIRCLE APPEARS ON THE ORBIT EVERY 6,2»j RADIANS OF ",/, 
"SOLUTION AND IS EQUIVALENT TO A KEY PHASOR MARK ON A CRO TRACE",/# 
"THIS IS TRUE ONLY WHEN TMAX IS A MULTIPLE OF 6,28,",/, 

"A PLUS SIGH APPEARS AT THE POINT ThE SOLUTION IS I N I T JALLY STARTED 
"» /» "OR CONTINUE) WITH RPCS > 1, "#//, 

"WHEN N * 12, THE ABSOLUTE JOURNAL MOTION APPEARS AS A DASHED",/, 
"LINE AND THE RELATIVE M 0 TI0N APPEARS AS A SOLID LINE. ",//, 

"THE CROSS COUPLING TERMS FOR THE BEARINGS ARE NUT PRINTED OUT*,/. 
"ON THE PLOTTER OUTPUT 8UT THEY DO APPEAR ON THE LP OUTPUT, ",/, 

"AS A SUGGESTION YOU COULO PUT A NEGATIVE CASE NuMBER WHEN AND ",/, 
"IF THE CROSS COUPLING TERMS FOR THE BEARINGS ARt NOT ZERO ,",/, 
"THIS WOULD INDICATE TO LOOK AT THE LP 1-0 FOR THE VALUES">)| 

ACAND I 

FMSHH * 0.0 J FMBHH ♦ 0.01 

yv ♦ i j 

RPC ♦ 1! 

READ (CR,/,TMAX,H,N,CASNO)[ALLD0NE)J 
READ (CR#/#RPM,N,K,OS,C I #QAC#E U#WJ, W 1 )) 
READ(CR,/#KXX#kYY,CXX,CYY,KXY»KYX»CXY#CYX! J 
IF N > 8 THEN READ(CR,/."K1X,KIY,CIX,C1Y)I 
READ(CR,/#FOR 1*0 STEP 1 UNTIL N DO (YtO#!!!)! 

WRITE (LPIPAGE I ) J 
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r# 

1 




NRITE(LP»<"INITIAL CONDITIONS"*)! 

NRITECLP,<F11.6>,FQR I ♦ 0 STEP 1 UNTIL N DOtYCO*!!))! 

REAO (CR,/,C$,ISCALE,XMIN,0X,XMIN2,DX2,XM!N3,DX3)1 

ir CS M THEN 

REAO(CR,/,RPCS)l 

RAO 4 RPM » 0,1047 I 

H *■ W / 386 J 

Ml 4 Ml/306 I 

MJ ♦ N J/ 386 ,01 

U 4 MX RADI 

U2 4 ZZxRAOf 

Z1 ♦ DS/ZZ1 Z24CI/ZZI Z34X/ZZ2I Z440AC/ZZ2! 

ZS, M/M J J Z6*(Cl42*CXX)/(?xZZ)7 ZZ24Z2/2! 

Z74(?xKXX4K)/(2xZZ2)| Z324Z3/2I 
Z84(Cl4 2 xCYY)/(2xZZ)| Z94(2XKYV4K)/(2XZZ2)J 
IF N > 8 THEN 

Z104 M/Ml I Z114K1X/ZZ?! Z124C1X/ZZI Z134KXX/ZZ2! 

Z144CXX/ZZI Z154K1Y/ZZ2I Z164CIY/ZZI 
Z174KYY/ZZ2I Z184CYY/ZZ7 

Z1S4KXY/ZZ21 Z204CXY/ZZI 7214KYX/ZZ2J Z234CYX/ZZI 

K?X 4 C2xKXXxkx(K42xkXX)44xKxCRA0xCXX)*2)/((K42xkXX)*2 ♦ 

«x(RA0xCXX)4 2 )| 

K2Y 4 (2xKYYxKx(K*2xKYY) 4 4xKx(RA0xCYY)*2>/ <(K42xKYY>*2 4 
4x(RADxCYY)*2)| 

C2X 4 2xKxk*CXX/((2xKXX 4K)*2 4 Ax (R AO“CXX 1*2) + DS 1 

C?Y 4 2XKXKXC YY/< (2xKVY*K)» 2 4 4x( RAO*C YY)*2 ) 4 OS I 

WCX 4 S«RT(2xKxKXX/((2xKXX4K)xM)> 1 

WCY 4 SgRT(2xKxKYY/((?XKYY4K)xM)) 1 

AC X 4 K2X/(WCXxC2X) l 

ACY 4 K2Y/(NCYxC2Y) | 

WCX 4 WCX/O, 10471 HCY 4 MCY/O, 10471 

TMAXH 4 TMAX 1 

Tm A X 4 TMAX 4 Y 1 0 ,0 ] 1 

REpfATl RPC 4 RPC 4 11 

TlME$ Tr P (TMAX, H» N, AY, NA, Y ) 1 

WRITE CLPtPAGE])! 

WRITE (LP,<X35, ’•ORBITAL POTION OF THE SINGLE MASS UNBALANCED ", 
"KOTOR", ///>)! 

WRITEUP, <X2, "CASE NO. ",F11*4.»//,X2, "ROTOR WEIGHT • *, 

F9,3," LB, ",X12, "ROTOR SPEEO «",F10,2»" RPM",/,X2, 

"ROTOR STIFFNESS «",F11,4," LB/MIL.",X3, 

"UNBALANCE »*,F8,J," MILS",/,X2, 

"SHAFT DAMPING *",F0,3»* LB-SEC/IN", X6, 

"INTERNAL DAMPING "",F9,3," LB-SEC/IN", /,X2, 

"CROSS COUPLING «",F9.2, " LB/IN", //,X2, 

"BEARING ST I F FN ESS", X 11, "BEARING DAMP ING",/,X3, 

"W XX »",F10,3," LB/MIL", X6,"CXX "",F10,3," LB"SEC/ IN",/, X3, 

"KYY *"»F 10,3," LB/MIL", X6,"CYY «",F10.3»" LB"SEC/IN">, 

C ASND,W, RPM# K/ 1000,0, EU, OS, C I, «AC#KXX/ 1000, 0,CXX,XYY/ 1000,0/ CYY)1 
HRITE(LP,<X3,"KXY «",F10,3," LB/MIL", X0,"CXY x",F10,3," LB-SEC/IN", 
/,X3> "K YX »",F10,3," LB/MIL", X6,"CYX »",F10,3#" LB-SEC/1N"> 

,KXY/1000,0,CXY,KYX/1000,0,CVX)1 

MRlTf (LP#<X3, "WEIGHT OF EACH BEARING "",F9,3," LB.",//>,WJ)1 
IF N > 0 THEN 

WHITE aP,<X2, "SUPPORT STIFFNE SS", X 1 1, "SUPPORT DAMPING", /,X3, 

" K 1 X F10 , 3» " LB/MIL'*»X6,"C1X »",F10,3," LB"SEC/IN"#/»X3» 


I 
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"K1Y •”,F10,3»" L8/MIL"#X«#"C1Y «"#FiO,3#" LB"SEC/IN"#//># 
K!X/t000,0#ClX#KlY/i000«0»ClY)3 

IF N > 6 THEN HR ITE(IP#<X3#"WEIGHT OF EACH SUPPORT «"#F9,3# 

• LB."#//>#Wn3 

WRITE(LP#<XA#"THAXH«"#F8,?#/,X7#"N • 13# /#X7#"H >*»rt,3*/«X4< 

"RPCS •"#IA>#THAXH#N#H#RPCS)3 


IF CS * 0 THEN 

begin 

IF VV • 1 then BEGIN 

PLOT (2#0#»A)3 PLOT <2#0#«5)3 


FILL 

ALPHAim NITM "Y-DIR,"#" 

tHIL"# 

"S3 "3 

FILL 

ALPHA2IO HITH "X-OIR,"#" 

CH IL"» 

"S3 "3 

fill 

alpi 1*3 

WITH 

"ABSOLU"#"TE 

ROT"#" 

OR MOT »# "ION "| 

FILL 

ALP2 t*J 

WITH 

"JOURNI 


MOTI*#" 

ON "I 

FILL 

ALPS I*) 

WITH 

"SUPPOR 


MOT I"#" 

ON "I 

FILL 

ALPA t*3 

WITH 

"N ■ 

f* # ft 

•#" 

RPM "1 

FILL 

ALP3 t* 3 

WITH 

’’W2 * 

it 

"#" 

LB, "3 

FILL 

ALP6 1*3 

WITH 

"EU ■ 


"#" 

M 1 1 * , "S "| 

FILL 

ALPF [*I 

WITH 

"KS • 

ft, ft 

"#" 

LB/*# l ’MIL ", 

FILL 

ALPB t *) 

WITH 

"WJ ■ 

f» # ft 

"#" 

LB, "1 

FILL 

ALP9 C *3 

WITH 

"CS ■ 

ft |ti 

"#" 

LB-S"»"EC/IN "3 

FILL 

ALPlOt*3 

WITH 

"KXX ■ 

ff # ft 

"#" 

LB"# "/M IL "3 

FILL 

ALPI 1 1 *3 

WITH 

"Cl » 

ff # ft 

"#" 

LB«S"#"EC/IN "| 

FILL 

ALP 121 *3 

WITH 

"KYY • 

tt # ft 

"#" 

L8"#"/MIL "3 

FILL 

ALPISI *3 

WITH 

m Q m 

ft jtt 

"#" 

LB/IN"3 

FILL 

ALP 1 4 E *3 

WITH 

"cxx - 

m p ft 

"#" 

LB*"#"SEC/IN"3 

FILL 

ALP 15 1 * 3 

WITH 

"WCX * 

ft # ft 

"#" 

RP"#"M "j 

FILL 

ALP16I*3 

WITH 

"CYY ■ 

ft jtt 

"#" 

LB-"#"SEC/IN"3 

FILL 

ALP 1 r t * 3 

WITH 

"WCY • 

ft, ft 

"#" 

RP"#"M "3 

FILL 

ALPI 6 C * 3 

WITH 

"ACX ■ 




FILL 

ALP20t*3 

WITH 

"ACY ■ 




IF N 

> 6 THEN 

BEGIN 




FILL 

ALP19 [ * 3 

WITH 

"W! • 

ft ,ft 

m f m 

LB "3 

FILL 

ALP21I*] 

WITH 

"K1X « 

ft # ft 

n t n 

L8/"#"MIL "3 

fill 

ALP22I *3 

WITH 

•cix • 

ft, ft 

m t n 

LB-"#"SEC/IN"3 

FILL 

ALP23I*] 

WITH 

"K1Y ■ 

ft , ft 


LB/"#"MIL "3 

FILL 

ALP2«t*3 

WITH 

"Cl Y • 

ft, ft 

n f * 

LB»"#"SEC/IN"j 

END J 







FILL 

ALP25t*3 

WITH 

"CYCLES 

tt|H 

m t m 

T"#*HR0UGH"3 

FILL 

ALP26 [ * 3 

WITH 

"FU • 

* p ft 

m f m 

LB,"#" "3 

FILL 

ALP2TI*] 

WITH 

"TRDB ■ 

ft, ft 

m 

»" ( "#" 

FILL 

ALP28E *3 

WITH 

"TRDS «"#" 

m 

#" ( "#" 

FILL 

ALP29[*3 

WITH 

"FB » 

•»,»» 

m f m 

LB, "3 

FILL 

ALP30I *] 

WITH 

•NO, 

"i 




END I 

W * VV ♦ 1 | 

BEGIN 

REAL KIRKRG I 

FOR I ♦ 0 STEP 1 UNTIL N DO 
YtO#I 3 ♦ AY [ I »NA ] 3 

NNITECLP»<X 2 #"END CONDITIONS FOR RPC ■ "# I 3 >#RPC* 1>3 

WRITE<IP#<FU, 6 >#FQR ! ♦ 0 STEP I UNTIL N Q 0 tT(O#I 3 I )3 

IM > 8 THEN BEGIN 

FOR ; ♦ 1 STEP 1 UNTIL NA 00 BE 6 IN 

BXRIP «. AY[5»l3*AYt9#U3 


)"3 

)*3 
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BYRtU ♦ AY£7#llfAYlllfI)J 

ENOi 

ENDJ 

IF I SCALE THEN BEGIN 
IF VV S 2 THEN BEGIN 
XMAX ♦ YMAX ♦ 3 x OX I 
YMIN ♦ XMIN f OY ♦ OX f 

XMAX? ♦ YHAX2 ♦ 3X0X21 
XMAX3 ♦ YMAX3 ♦ 3*0^31 
YMIN? * XMIN2 I DY2 ♦ DX2| 

YMIN3 ♦ XMIN3J DY3 «• 0X3| 

ENDI 

PL0TCHEK(AY[1»*J/ AY13#*1i»NA/«XNIN,XMAX»YHIN,YHAX)| 

SCALES (AY(3f*]*NA#YMlN»0Y»CS)J 
SCALES !AYil,*)fNA#XMIN#DXfCS)l 

PL&TCHEK(AYI5»*J#AYt7»* J#NA»XMIN2»XMAX2#YMIN2»YMAX2)I 
SCALES! AY! 5f* If NAfXMIN2,DX2fCS>| 

SCALES! AYl7f *)#NAf YMlN2f DY2fCS)| 

IF N > 8 THEN CEGIN 
I ♦ 1U 

PL0TCHEK(AYC9#*]/AYtI»*]#NA#XNlN3#XMAX3»YMIN3» YMAX3JI 

SCALES (AYI I #*1#NA»YMIN3#0Y3#CS)I 

SCALES(AYt9#*l»NA*XMIN3#DX3*CS)| 

PL0TCHEK(BXR*BYR,NA,XMIN2,XMAX2#YMIN2»YMAX2)J 

SCALES(BYR#NA»YMIN2»DV2#CS)J 

SCALE S(8XR#NA#XMIN2,DX2#CS)I 

ENDI 

END 

ELSE 

BEGIN 

SAMESCALE(AYI1#*J»AY13»*]#NA#XMIN#XMAX,DX»YMIN#YMAX»DY)I 
PLOTCHEK(AYCl,*],AYt 3,*]»NA#XMIN,XMAX#YMIN,YHAX)| 

SAME SCALE! A Y[ 5f * 1 f AY[ 7, * ]f NAf XMIN2f XMAX2f DX2f YMlN?f YMAX2f 0Y2)> 

PLOTCHEK! AYI5»*J# AY[7»*J#NA#XMIN2»XMAX2#YMIN2»YNAX2)I 

SCALES ! AY!3f *JfNAfYMlN#DVfCS)l 

SCALES ! AY[lf *JfNAfXMlNfDXf CS)J 

SCALES! AYI7#*J#NA»YMIN2#i)Y2#CS) f 

SCaLES!AY15>*]»NA*XMIN2#0X2#CS)I 

IF N > 8 THEN BEGIN 

PL0TCHEK!PXH*BYR*NA,xMIN2»XMAX2»YMIN2«YM/4X2)| 
SCALES!BXR»NA#XNIN2#0X2#CS)I 
SCALE S!BYR*NA#YMIN2>DY2#CS)I 
1 ♦ 111 

SAMESCaLE<AYC9 # *),AY[1,*J,NA,XMIN3,XMAX3#DX3,YMIN3,YNAX3,DY3)J 

PL0TCHEK(AY(9»*]#AYtI>*]fNA>XMIN3»XMAX3»YNIN3*YNAX3)l 
SCALES! AY !I/*1»NA»YHIN3»DY3#CS)I 
SCALES! AYt 9#* 1#NA*XNIN3#0X3#CS)| 

ENDI 
END | 

IF VV S 2 THEN 
BEG T N 

OMBITTOPIDI 

AGRID!XMIN#DX#YMIN#DY1) 

END J 

SYMBOL !AY[ 1# 11* AY[3»ll« ( lA« ALFHA1»0*-1 3)1 
XA ♦ 6,26/HI 


i 
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FOR I * X A STEP XA UNTIL NA DO 
SYNB0L(AYtl#I)#AYtJ#lJ#0,07#ALPMAl#0#"5)l 
LYNEt AYt l»*).'AYt3»*)#NA#CS)) 

PL0T(12»0#"5)I 
If VV S 2 THEN 

begin 

ORBlTTO.*'<?)J 

A«RI0(XNIN2/0X2»YNIN2»0Y2)| 

end; 

$VHB0L(AYt5#n»AYt7#i)#0'!4»ALPHA!»0»M3)f 

FOR I ♦ XA STEP XA UNTIL NA DO 

SYNBOLt AYt5# I )#A Yt 7» I )#0»0F# ALPHA1 »0#*5 )l 

LYNE(AY(5»*J#AYtr»*J»NA#CS)l 

If N > 0 THEN BEGIN 

SYNBOKBXRIIJfBYRII J,0,U# ALPHA! #0#»1 3)1 

FOR I * XA STEP XA UnTI. N A DO 

SYMBOL t§XRU)#BVRtl)#0. 07# ALPHA1,0#-5)I 

DASHl INE(8XR#BYR»NA#CS)I 

ENpj 

PLOT ( 12»0»»5) > 

I r N > 0 THEN BEGIN 
If VV S 2 THEN 
BEGIN 

ORBITTOP(S)! 

AGRID(XNIN3#DX3#YNIN3»0Y3); 

END) 

J ♦ 111 

SYMBOL (AYt9# 1 1# AY tU#t)#0«lG# ALPHA! #0#"1 3)1 

FOR I + XA STEP XA UNTIL NA DO 

SYNBOL ( AYI9# U# AYN# IJ#0, 07# ALPMAl>0#-5) | 

LYNEf AYI9#* J# AYt J»* 1#NA#CS)I 
ENOl 

ENOl 

NRIU(LPCOBU)l 

L ITFCLP.-<"TRDB-HAX » "»F11.4»" AND OCCURS AT "# FT, 2# 

» Cf CLES*>» IF EU»0.0 THEN S0RTCFNBH)>«2,0 EL?E SORT (FNBH)*2.0/EU» 
TIHHB/0,2 0)1 
IF N>0 THEN 

WRITE(LP#« ,, TRDS-HAX « *#F! 1 »* ■ "AND OCCURS AT "#F7,2# 

" CYCLES"*# IF EU«0« 0 THEN SQRT(FMSH)*2.0 ELSE S0RT<FMSH)K2*0/EU# 
TINHS/0.20)! 

NRITECLP#<"UNBALANCE STIFFNESS • *#Flt,2# " LB/MJL"#X2# 

"AND UNBALANCE • "#F0»3#" N1LS"*#ZZ2/I000,0#EU)J 
IF FMBH > FHBHH THEN BEGIN 
fnbhh *■ FMBHJ tinmbh ♦ TIMHBI 

END) 

IF FMSH > FMSHH THEN BEGIN 
FNSHH «■ FHSHI TINHSH ♦ TINHSI 
END| 

IF RPC S RPCS THEN 
BEGIN 

PL0T("24#O#"5)J 
TMAX ♦ TMAX ♦ TMAXH) 

ISCALE ♦ TRUE I 
GO TO REPEAT) 
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I 


ENp| 

IF EU • 0.0 THEN BEGIN 
FM6HH ♦ S8RT!FMBHH)*2.0I 
FMSHH «• S0RT(FM$HH)k2,0I 
END ELSE 
BEGIN 

FMBHH «■ S0RT(FMBHH)*2,0/EUI 
FMSHH ♦ SfiRT(FNSHH)K2*0/EU| 

END) 

TlMHBH ♦ T I HHBH/6 *26 J 
TINHSH ♦ TIMHSH/6.281 
PlOT!*24,0,-5)) 

IF N>6 THEN BEGIN 

SYMBOL !0 *75# 6 ,50,0.14# ALP29#0# 16)1 
NUMBER! 1. 23.6 . 50,0, 14 ,FMBHHxZZ2xEU/ 1000,0 ,0,3)1 
ENOI 

NUMBER! 4, 31 #6.75# .I4,AY(0,NA)/6.28,0#2)) 

SYMBOL (3 .75 #6, 50#. 14, AtP26*0#21 ) J 
NUMBER (4,23# 6. 50, .14,722/ 1000, 0*EU,0# 3)) 

PLOT ! 12#0#**5 ) I 
IF N*6 THEN BEGIN 

SVM80L(0,75*6.50,0,14»ALP29#0#16)I 
NUMBER! 1.23*6.50,0. 14»FMBHH*ZZ?*EU/ 1000, 0* 0, 3) I 
END) 

NUMBER(4.31*8.75#,14,AYt0#NA)/6,26,0,2)) 

SYMB0L(3,75#6,50# ,14, ALP27, 0,24)) NUMBER (4, 35, 6, 50#, 14# FMBHH# 0,31) 
N"MBER(5.43,6.50,,14,TIMHPH#0,2)) 

PLOT ( 12#0#*5)I 
IF N>8 THEN BEGIN. 

NUMBER(4.3l#6.75,*14,AYtO#NA}/«. 28,0, 2)1 

SYMBOL (3, 75, 6. 50, ,14, ALP26#0# 24 )|NI'M6ER!4. 35,6,50 *,14, FMSHH, C, 3)1 
NUMBER(5.43,6.50#.14,TIMHSH#0#2)) 

END) 

PLOT! 12,0, *5)1 
GO TO ACARD ) 

ALLOONE I 
PLOT ! 1# C -3)) 

ENO) 

WRITE !LPtPAGE))l 

WHITE !LP»<"T0TAL PROCESSOR TIME « "# F6,2#*1#"MINUTES">, 

TIME!2) / 3600 ) ) 

WHITE ILPtPAGE ) , <*TOTAL 1*0 TIME ■ "# F6,2* XI, •MINUTES* > # 

T I ME! 3 ) / 3600 ) I 

ENO. 

ARCTAF IS SEGMENT NUMBER 0080, ?RT ADORESS 1$ 0117 

COS IS SEGMENT NUMBER 0061, PRT AODRESS IS 0075 

EXP IS SEGMENT NUMBER 0082, PRT AOORESS IS 0072 

LN IS SEGMENT NUMBER 0083*PRT AOORESS IS 0071 

SIN IS SEGMENT NUMBER 0084, PRT AOORESS IS 0076 

SORT IS SEGMENT NUMBER 0065, PRT AOORESS IS 0436 

OUTPUT !W) IS SEGMENT NUMBER 0086#PST ADDRESS IS 0440 

BLOCK CONTROL IS SEGMENT NUMBER 00B7#PRT AOORESS IS 0005 

INPUT ! W ) IS SEGMENT NUMBER 0066«PRT ADDRESS IS 04AA 

X TO THE I IS SEGMENT NUMBER 0089, PRT ADDRESS IS 00T3 

GO TO SOLVCR IS SEGMENT NUMBER 0090#PRT AOORESS IS 0065 

Algol write is segment number oo9i,prt aooress is ooia 
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ALGOL REAP IS SEGMENT NUMBER 0092#PRT ADDRESS IS 0015 

ALGOL SELECT IS SEGMENT NUMBER 009S»PRT ADDRESS IS 0016 

FILE ATTRBUTS IS SEGMENT NUMBER 009A.PRT ADDRESS IS 0033 

COMPILATION TIME(SECQNOS)! PR * 57 1/0 • 162 

NUMBER OP ERRORS DETECTED ■ 000* LAST ERROR ON CARO # 

NUMBER OF SEQUENCE ERRORS COUNTED ■ 0. 

NUMBER OF SLOM MARNINGS * 0, 

PRT SIZE* 33AJ TOTAL SEGMENT SIZE* 3952 NOROS. 

DISK STORAGE ReQ«* 360 SeGS.) NO, SEGSt* 95, 

ESTIMATED CORE STORAGE REQUIREMENT • 15355 NOROS, 
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PLOTTER OUTPUT INFORMATION AND SUGGESTIONS 

MAKE TMAX A MULTIPLE OF 6.28 (BUT LESS THAN 35 IF HaO.OS) 

A REASONABLE VALUE OF H IS 0*05 WHICH 6IVES 125 STEPS PER 
CYCLE OF RUFNING SPEED FOR THE INTEGRATION. 

TWO C2) MINUTES PROCESSOR TIME IS REQUIRED FOR 10 CYCLES OF 
SOLUTION FOR N • 8. 

THREE (3) MINUTES PROCESSOR TIME IS REQUIRED FOR 10 CYCLES i)F 
SOLUTION FOR N • 12. 

A SMALL CIRCLE APPEARS ON THE ORBIT EVERY 6.26 RADIANS OF 

solution and is eouivaicnt to a key phasor mark on a C Rc trace 

THIS IS TRUE ONLY WHEN TMAX IS A MULTIPLE OF 6.26. 

A PLUS SIGN APPEARS AT THE POINT T«E SOLUTION IS INITIALLY STARTEO 
OR CONTINUED WITH RPCS > 1. 

WHEN N « 12# THE ABSOLUTE JOURNAL MOTION APPEARS AS A DASHED 
LINE ANO THE RELATIVE MOTION APPEARS AS A SOLID LINE. 

THE CROSS COUPLING TERMS FOR THE BEARINGS ARE NOT PRINTED OUT 
ON THE PLOTTER OUTPUT BUT THEY DO APPEAR ON THE LP OUTPUT. 

AS A SUGGESTION YOU COUlD PUT A NEGATIVE CASE NUMBER WHEN AnD 
IF l HE CROSS COUPLING TERhS FOR THE BEARINGS ARE NOT ZERO . 

THIS NOULO INDICATE TO LOOK AT THE LP 1-0 FOR THE VALUES 
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ORBITAL MOTION OF TMC SINGLE NASS UNBALANCED MOTOR 


CASE NO, 707,7100 

ROTOR HEIGHT ■ 675*000 LB* 

ROTOR STIFFNESS " 260,0000 LB/NIL* 

SHAFT OAMPING • 6.000 L9*S£C/lN 

CROSS COUPLING ■ 0,00 LB/IN 

bearing stiffness 

KXX • 351*000 LB/MlL 

KYY * 606*000 LB/MIL 

KXY * 0*000 LB/MIL 

KYX • 0.000 LB/NU 

HEIGHT OF EACH BEARING • 


TMAXH« 12,56 
N • 8 

H ■ 0,050 

RPCS * 5 

ENO CONOITIONS FOR RPC » 1 

75*400000 
•0,756303 
•0.497459 
1,746020 
•0. '37636 
0,366998 
-0*206504 
0,435169 
0,060061 


ROTOR SPCEO • 10000,00 RPM 

UNBALANCE • 0,500 NILS 

INTERNAL OAMPING • 100*000 LB-SEC/IN 


BEARING DAMPING 
CXX ■ 739,000 LB-SEC/IN 

CYY « 665.000 LB-SEC/IN 

CXY • 0.000 LB-SEC/IN 

CYX ■ 0.000 LB-SEC/IN 

312,000 LB, 


TRDB-MAX ■ 0,7590 ANp OCCURS AT 11,82 CYCLES 

UNBALANCE STIfFNESS « 1916.95 LB/NIL AND UNBALANCE ■ 0.500 MILS 


INITIAL CONDITIONS 
62.600000 
•1,690000 
0,163960 
•0,601630 
•1*133000 
•0,353260 
>0,167480 
0,126290 
-0.107260 
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NOMENCLATURE 

Amplification factor at rigid support critical 
Complex bearing amplitude, in 
Complex support amplitude, in 
Complex rotor amplitude, in 
Damping ratio = C 1 /C 2 (DIM) 

Bearing damping, Ib-sec/in 

Critical damping coefficient, Ib-sec/in 

Rotor internal damping, Ib-sec/in 

Absolute shaft damping, Ib-sec/in 

Support damping, ib-sec/in 

Effective rotor-bearing damping, Ib-sec/in 

Rotor mass eccentricity, in 

Force transmitted to foundation, lb 

F orce transmitted to bearing housing, lb 

Stiffness ratio, Ki/K 2 

Bearing stiffness, lb/ in 

Rotor-shaft stiffness, lb/ in 

Support stiffness, lb/ in 

Effective rotor-bearing stiffness, Ib/in 

Mass ratio, = M!/M 2 (DIM) 

Support mass, lb-sec 2 / In 

Rotor mass, lb-sec 2 / in 
Rotor critical speed, CRFW!] 

1st node point on response plot 
Rotor cross-coupling stiffness, Ib/in 
2nd node on response plots 
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V * 2 


(DIM) 
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Rotor absolute displacement amplitude 
Kinetic energy 

Transmissibi I it/ = Fi/(M 2 e u « 2 ) 

Potential Energy 
Velocity, in/sec 
Support weight, lb 
Defined a?. X 2 

Shaft re.ative displacement in x - direction 
Support displacement in x - direction 
Rotor absolute displacement in x - direction 
Journal relative displacement in x - direction 
Shaft relative displacement in y - direction 
Support d i sp I acement i n y - d i rect i on 
Rotor absolute displacement in y - direction 
Journal relative displacement in y - direction 
Complex shaft relative amplitude, in. 

Complex support amplitude, in. 

Comp I ex rotor amp I i t ude , in. 

Complex journal amplitude, in. 

Rotor angular acceleration, rad/sec 2 

Phase angle of support motion relative to rotor unbalance, DEG 
Phase angle of rotor motion relative to rotor unbalance, DEG 
Phase angle of bearing motion relative to rotor unba:ance, DEG 
Angular displacement, rad 
Defined as K/M (DIM) 

Damping ration= Ci/c c (DIM) 

Rotor absolute amplitude phase angle, deg. 


T 
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♦ 

X 

? 

Ill 

“ 1,2 

u 

c 

ill, fl 2 


A 

CB 

CD 

DC 

E 

FTRSWC 

N 

QAC 

TRDB 

TRDS 

FU* 

W 

WC 



Moment of inertia 

Optimum amplitude for tuned system 

Defined as Rj 2 or 0 2 2 w ^en calculating required damping at point 
P or Q respectively 

Ftotor angular velocity, rad/ sec 

Rotor system critical speeds, rad/sec 

Rigid support critical speed, rad/sec 



Speeds at which the node point P and 0 occur on response plots 


(AUTOMATIC PLOTTER NOMENCLATURE) 

Amplification factor at rigid support critical (DIM) 

Bearing damping, lb-sec/ in 

Shaft damping coefficient, lb-sec/ in 

Internal damping, lb-sec/ in 

Rotor mass eccentricity, in 

Force transmitted at rigid support critical 

Rotor speed, RPM 

Aerodynamic cross-coupling coef-, Ib/in. 

Maximum bearing force transmitted (DIM) 

Maximum support force transmitted (DIM) 

Rotating unbalance ioad per mil unbalance eccentricity, lb 

Rotor speed, rad/sec 

Rigid support critical speed, rad/ sec. 
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CRITICAL FREQUENCY VS SPRING RATE 
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BEARING SUPPORT STIFFNESS KATIO, K=K,/K 2 (DIM.) 

Figure 2 Dimensionless Critical Speeds vs, Suppor+ Stiffness 
Ratio for Va rious Sup port Hous ing Mass Ratios 
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ROTOR AMPLITUDE AT CRITICAL vs. MASS RATIO 
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Figure 13 Ro+or Amplitude at Critical Speeds vs. Mass Ratio 
Various Values for Various Damping Ratios, A = 10 








ROTOR AMPLITUDE AT CRITICAL vs. MASS RATIO 



Figure 15 Rotor Amplitude at Critical Speeds vs. Mass Ratio 
for Various Damping Ratios, A = 10, K = 5.0 










ROTOR MAXIMUM AMPLITUDE VS DAMPING RATIO 
FOR VARIOUS STIFFNESS RATIOS 
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Figure 16 Rotor Maximum Amplitude vs. Damping Ratio for 

Various Values of Stiffness Ratios for a Low Mass 
Ratio Support, M = 0.01, A = 10 










Figure 18 Rotor Mrx I mum Amplitude vs. Dampl 





Figure 19 Rotor Maximum Amplitude Vs Damping Ratio for Various Values of Stiffness 
Ratios for a High Mass Ratio Support, M = 2, A « 10 




ROTOR MAXIMUM AMPLITUDE FOR VARIOUS VALUES OF STIFFNESS AND MASS RATIO 

WITH OPTIMUM SUPPORT DAMPING 



Figure 20 Rotor Maximum Amplitude for Various Values of Stiff- 
ness and Mass Ratio with Optimum Support Damping 






ROTOR CRITICAL AMPLITUDE vs. SUPPORT DAMFING RATIO FOR VARIOUS VALUES OF 

SUPPORT STIFFNESS WITH OPTIMUM MASS RATIO 
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C.4 0.6 0.8 I 2 4 

DAMPING RATIO, C*C,/C 2 (DIM) 

Figure 22 Ro+or Critical Amplitude vs. Support Damping Ratio 
for Various Values of Support Stiffness with Optimum 
Mass Ratio 
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Figure 23 Optimum Damping and Mass Ratios for Various Values 
of Stiffness Ratio 
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RBSCLUTE ROTOR MOTION 


N - 30000 RPR 

H ■ 

0. 100 

K - 0.092 

C • 

43.946 

HZ « 96.60 L6. 

KB * 

SX,0X LB/IN 

KS - 500,000 IS/ IN 

GB « 

100.0 LB-SEC/IN 

X • 0.S LB-SEC/IN 

HI - 

9.66 LB. 

X » 0.S LB-SEC/IN 

K1 * 

25,000 LB/ IN 

(MC - 0.5 IB/IN 

Cl * 

1000.0 LB-SEC/IN 

TUB - 0.585 AND OCCURS AT 

0.58 CYCLES 

TRDS « 2.159 RNO OCCURS RT 

0.99 CYCLES 

FU«- 2469.026 LBS 

• 



UBO — *4. ODD 


- 2.000 0.000 
X-DIR. 


2.000 


Figure 24 Dimensionless Transient Motion of an Unbalanced Rotor 
for Twelve Cycles on Over-Damped Supports C« = M * 
0.1, C = 44j 
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BEARING MOTION 

N ■ 30000 RPH M - 0.100 

K - 0.092 C - 93.946 

M2 * 96.60 LB. KB - 500,000 LB/IN 

KS « 500,000 LB/ IN CB - 100.0 LB-SEC/ IN 

DC • 0.5 LB-SEC/IN Ml • 9.66 LB. 

CO - 0.6 LB-SEC/IN K1 - 25,000 LB/IN 

QRC « 0.5 LB/IN Cl = 1000.0 LB-SEC/IN 

TR08 « 0.585 AND OCCURS RT 0.58 CYCLES 

TRDS ■ 2. 159 RNO OCCURS RT 0.99 CYCLES 



•- 3.000 - 2.000 - 1.000 0.000 1.000 2.000 3.000 


X-UIfi. 

Figure 25 Dimensionless Bearing Absolute and Relative Transient 
Motion for Twelve Cycles on Over-Damped Supports 
[K = M = 0. I , C * 44] 
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SUPPORT MOTION 




L 3.000 - 2.000 - 1.000 0.000 1.000 2.00 0 3.000 

X-DIR. 


Figure 26 Dimensionless Transient Support Motion for Twelve 
Cycles wlfh Excessive Damping [K * M * 0.1, C * 44} 
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ABSOLUTE ROTOR MOTION 


N - 

30000 RPM 

H « 

0.100 

K « 

0.092 

C » 

0.H39 

M2 « 

96.60 LB. 

KB « 

500,000 LB/IN 

KS - 

500,000 LB/IN 

CB - 

100.0 LB-SEC/ IN 

DC * 

0.5 LB- SEC/IN 

Ml - 

9.66 LB. 

CD - 

0.5 LB-SEC/ IN 

K1 « 

25,000 LB/IN 

QRC - 

• 99.0 LB/IN 

Cl - 

10.0 LB-SEC/ IN 

TRDB 

- 0. 16H RND OCCURS RT 

8.81 CYCLES 

TRDS 

* 0.088 RND OCCU 

RO RT 

8.09 CYCLES 


2469.028J3S 



000 - B . 0 C 0 -* 1.000 0.000 * 1.000 0.000 12.000 

X-DIR. 

Figure 27 Dimensionless Transient Motion with Under-Damped 
Flexible Supports for Twelve Cycles CK * M = 0.10, 

C = 0.443 
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Y-DIR 
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BEARING MOTION 

N - 30000 RPK 

K • 0.092 
H2 - 96.60 LB. 

KS « 500,000 LB/IN 

OC - 0.5 LB-SEC/IN 

CD - 0.5 LB-SEC/IN 

QflC - 99.0 LB/IN 





H ■ 0.100 

C - 0.439 

KB * 500,000 LB/IN 

CB - 100.0 LB-SEC/IN 
HI - 9.66 LB. 

K1 » 25,000 LB/IN 

Cl • 10.0 LB-SEC/ IN 


Figure 28 


Dimensionless Bearing Absolute and Relative Transient 
Motion for Twelve Cycles on Under-Demped Supports 
[K - M ■ 0.1, C « 0.44] 
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SUPPORT MOTION 



“i ~i r 

^ 2.000 - 8.000 - 4.000 0.000 

X-DIR. 


4.000 


6. 000 12.000 


N • 

30000 RPM 

H « 

0.100 

K « 

0.092 

C * 

0.439 

HZ - 

96.60 LB. 

KB « 

500.000 LB/IN 

KS » 

500.000 LB/IN 

CB - 

100.0 LB-SEC/IN 

DC - 

0.5 LB-SEC/IN 

HI « 

9.66 LB. 

CD * 

0.5 LB-SEC/ IN 

K1 * 

25.000 LB/IN ! 

QRC • 

* 99.0 LB/IN 

Cl - 

10.0 LB-SEQ/IN ! 

TRDB 

« 0. 164 AND OCCURS AT 

8.81 CYCLES 

TRDS 

- 0.088 (WO OCCURS RT 

8.09 CYCLES 1 


FU*« 2469. 026LBS 

i 

t — 


Figure 29 Dimensionless Transient Support Motion for Twelve 
Cycles with Under Damping [K * M = 0.1, C ® 0.44] 
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ABSOLUTE ROTOR MOTION 


N - 30000 RPM H « 

K - 0.092 C - 

H2 - 96.60 LB. KB 

KS • 500,000 LB/IN CB 

DC • 0.5 LB-SEC/IN HI 

CO • 0.5 HVfBC/IN K1 

one • o.«^wir ci 

TROB « 0.266 AC-O&URS RT 


0.100 

SOO.ioOO LB/IN 

100.0 LB-SEC/IN 
9.66 LB. 

2S,000 LB/ IN 

125.0 LB-SEC/IN 
0.19 CYCLES 



Figure 30 Dimensionless Rotor Motion with Optimum Steady-o . ate 
Damping Showing the Steady-State Orbit After Seven 
Cycles of Running Spepd f: , = M = 0.1, C = 5.5J 
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BEARING MOTION 

N - 30000 RPH H > 0.100 

K * 0.092 0 - 5.193 

M2 - 96.60 LB. KB • 500.000 LB/IN 

KS • 500.000 LB/!N CB « 100.0 LB-SEC/IK 

DC • 0.5 LB-SEC/IN Ml » 9.66 LB. 

CD < 0.5 LB-SEC/IN K1 « 25,000 LB/ IN 

QRC - 0.5 LB/IN Cl - 125.0 LB-SEC/IN 

TRDB = 0.2f c FIND OCCURS FIT 0.19 CYCLES 



X-DIR. 


Figure 31 Dimensionless fearing Absolute and Relative Transient 
Motion with Optimum Steady-State na ': 1 '>!ng fi 1 ' ! M s 0.1, 
C « 5.5] 
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